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Abstract 


EXPERIMENTAL  INVESTIGATION  OF  A  SUPERSONIC  TURBULENT  BOUNDARY 
LAYER  WITH  ADVERSE  PRESSURE  GRADIENT 


Laser  Doppler  Velocimetry  (LDV)  measurements  were  made  to  quantify  the  effect  of  adverse 
pressure  gradient  on  the  compressible  turbulent  flow  structure  in  a  Mach  2.9  boundary  layer 
(Re/m=  1.75  x  10^).  Measurements  included  profiles  of  2-D  mean  velocities,  turbulence 
intensities,  Reynolds  shear  stresses,  intermittency,  flatness  and  skewness.  In  addition,  mean 
strain  rates  were  also  measured.  The  boundary  layer  measurements  were  acquired  for  both  flat 
plate  and  compression  ramp  models.  LDV  measurements  were  made  at  two  locations  on  the 
compression  ramp  model  at  68  cm  and  71  cm  downstream  of  the  nozzle  throat.  At  these 
locations,  P  =  1.12  and  -0.94,  respectively.  Flow  visualization  was  accomplished  by  nanosecond 
shadowgraph  and  schlieren  photography.  Results  indicate  that  the  adverse  pressure  gradient 
increased  the  Reynolds  shear  stresses  by  190%  and  turbulence  intensities  by  24%  of  the  flat  plat 
values.  In  the  favorable  pressure  gradient  region,  these  quantities  were  decreased  by  52%  and 
7%  with  respect  to  the  adverse  pressure  gradient  values,  respectively. 
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Experimental  Investigation  of  a 
Supersonic  Turbulent  Boundary  Layer  with 
Adverse  Pressure  Gradient 


1.  Introduction 

1.1  Background 

The  study  of  aerodynamics  has  been  improved  by  the  numerical  capabilities  of  present 
day  computers  to  model,  or  simulate,  naturally  occuring  flow  structures.  Computational  fluid 
dynamics  (CFD)  creates  a  mathematical  model  that  approximates  these  flow  structures.  By  using 
numerical  methods,  many  flow  structures  may  be  evaluated  without  the  necessity  of  building  a 
model,  and  performing  wind  tunnel  experiments.  Turbulent  flow,  which  can  be  characterized  by 
the  unsteady  Navier-Stokes  equations,  is  usually  present  when  flows  have  high  Reynolds 
numbers.  Due  to  the  large  variation  in  the  spatial  and  temporal  length  scales  of  the  turbulent 
flow  structure,  closed  form  solutions  to  practical  high  Reynolds  number  problems  are  not 
possible.  Current  computational  capabilities  are  not  able  to  handle  these  problems  within  a 
reasonable  amount  of  time.  For  example,  flow  over  a  cylinder  of  Re^  =  lO’  would  require 

10^^  numerical  operations.  A  computer  that  is  capable  of  one  operation  per  nanosecond  would 
require  320,000  years  to  complete  this  analysis  [34].  Thus,  turbulent  flow  solutions  are  typically 
derived  from  time-averaged  approximations,  such  as  the  Reynolds  and  Favre  averaging  methods. 
The  time-averaged  form  of  the  governing  equations  is  solved  by  casting  the  dependent  variables 
into  the  sum  of  a  mean  and  fluctuating  component.  Because  of  the  non-linear  terms  in  the 
unsteady  Navier-Stokes  equations,  this  process  introduces  additional  cross-correlations  of 
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fluctuating  terms,  resulting  in  more  unknowns  than  governing  equations.  Turbulence  modeling 
attempts  to  relate  these  new  terms  to  other  flow  parameters  to  provide  enough  equations  to  close 
the  system. 

Wilcox  [35]  has  stated  that  “an  ideal  model  should  introduce  the  minimum  amount  of 
complexity  while  capturing  the  essence  of  the  relevant  physics.”  In  the  past,  researchers  have 
extended  empirically  based  incompressible  models  to  high-speed  flows  on  an  ad  hoc  basis.  This 
is  a  natural  progression  because  a  great  deal  of  incompressible  boundary  layer  research  has  been 
conducted.  The  use  of  these  incompressible  models  for  compressible  flow  fields  has  had  limited 
success  [3]. 

Turbulence  modeling  has  progressed  as  greater  knowledge  and  understanding  of 
turbulent  flow  structures  has  increased.  Turbulence  models  are  developed  by  making  both 
observations  and  assumptions  about  the  flow,  and  then  generating  a  model.  The  numerical 
results  are  then  compared  to  experimental  results  to  determine  the  validity  of  the  model.  If  the 
model  predicts  the  actual  flow  characteristics  quite  well,  then  it  is  extended  to  more  complex 
flow  situations  to  determine  its  robustness.  If  the  model  fails,  then  it  is  discarded,  and  further 
observations  and  assumptions  are  required. 

Thus,  accurate  experimental  data  are  required  in  order  to  generate  accurate  turbulence 
models.  There  has  been  extensive  research  into  the  nature  of  turbulent  boundary  layers,  both  in 
the  subsonic  and  supersonic  regimes.  Most  of  these  investigations  have  studied  the  turbulent 
boundary  layer  flowing  over  a  flat  plate.  A  great  deal  of  knowledge  has  been  gained,  but 
unfortunately,  aerodynamic  flows  encounter  more  than  just  flat  plates.  Of  particular  interest,  is 
the  study  of  turbulent  boundary  layers  with  pressure  gradient  effects.  The  effects  of  pressure 
gradients  tend  to  cause  large  changes  in  the  turbulent  flow  structure.  These  types  of  flow  are 
prevalent  in  both  internal  and  external  flows. 
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1.2  Objectives 

The  primary  objective  of  this  research  was  to  investigate  experimentally  the  importance 
of  adverse  pressure  gradients  (APG)  on  supersonic  turbulent  boundary  layers.  A  great  deal  of 
studies  published  in  the  literature  have  investigated  the  effects  of  adverse  pressure  gradients  on 
turbulent  boundary  layers  (see  Refs.  [4,  9,  13,  22,  24,  29,  33]).  These  studies  were  made  with 
conventional  probe  and  hot-wire  anemometry  techniques.  The  present  study  used  Laser  Doppler 
Velocimetry  (LDV)  to  measure  the  mean  and  fluctuating  velocities.  Several  other  researchers 
have  used  LDV  to  make  supersonic  turbulent  boundary  layer  measurements  over  a  flat  plate  (see 
Refs.  [12,  18-19,  25]).  After  an  extensive  literature  search,  it  was  determined  that  this  is  the  first 
investigation  to  use  modem  LDV  technique  as  a  means  to  investigate  a  supersonic  turbulent 
boundary  layer  subject  to  an  adverse  pressure  gradient.  This  clear  lack  of  data  was  also 
mentioned  by  Spina  [28]. 

A  second  objective  of  the  present  study  was  to  provide  accurate  empirical  data  to 
validate  the  turbulence  models  created  by  computational  fluid  dynamicists.  Accurate  data  are 
required  to  complete  this  task.  This  was  in  fact  accomplished,  as  Pick  [14]  has  validated  a  k-to 
turbulence  model  [35]  against  the  present  LDV  results. 

A  third  objective  is  to  compare  the  present  LDV  results  with  the  conventional  and  hot¬ 
wire  anemometry  results  generated  by  Dotter  [9]  in  the  same  AFIT  Mach  3.0  wind  tunnel  with 
the  same  compression  ramp  model.  In  this  manner,  the  validity  of  the  assumptions  used  to 
decompose  the  fluctuating  turbulent  terms  of  that  data  may  be  accessed. 

1.3  Settles  and  Dodson  Criteria 

This  experiment  was  conducted  with  the  intent  of  satisfying  the  criteria  established  by 
Settles  and  Dodson  [26]  to  provide  a  valid  database  of  experimental  results  for  turbulence 
modeling  and  CFD  code  validation.  This  criteria  is  summarized  below: 
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1 .  Baseline  Applicability:  All  candidate  studies  must  be  experiments  involving 
turbulent  flows  in  either  the  supersonic  or  hypersonic  Mach  number  range  (i.e.  M  ~ 
3.0  or  higher). 

2.  Simplicity:  All  candidate  studies  passing  this  criterion  must  involve  experimental 
geometries  sufficiently  simple  that  they  may  be  modeled  by  CFD  methods  without 
enormous  difficulty. 

3.  Specific  Applicability:  All  candidate  studies  passing  this  criterion  must  be  capable 
of  providing  some  useful  test  of  turbulence  modeling. 

4.  Well-defined  experimental  boundary  conditions:  All  incoming  conditions 
(especially  the  stat  of  the  incoming  boundary  layer)  must  be  carefully  documented. 
For  studies 

claiming  “two  dimensional”  flow,  data  indicating  the  extent  of  the  spanwise  flow 
variations  should  be  provided. 

5.  Well-defined  experimental  error  bounds:  The  experimenter  must  provide  an 
analysis  of  the  accuracy  and  repeatability  of  the  data,  or  error  bars  on  the  data 
themselves.  Further,  error  bounds  on  the  data  must  be  substantiated  in  a  quantifiable 
manner. 

6.  Adequate  documentation  of  data:  Data  must  be  documented  and  tabulated  in  a 
machine-readable  form. 

7.  Adequate  spatial  resolution  of  data:  Experiments  must  present  the  data  of 
sufficiently  high  resolution,  compared  with  the  scale  of  the  flow  in  question,  that  the 
key  features  of  the  flow  are  clearly  resolved. 


1.4  Synopsis  of  the  Present  Investigation 

Laser  Doppler  Velocimetry  (LDV)  was  used  to  make  mean  and  fluctuating  velocity 
measurements  in  the  AFIT  Mach  3.0  wind  tunnel.  The  measurements  were  made  over  flat  plate 
and  compression  ramp  models  at  the  same  distance  from  the  nozzle  throat.  In  this  manner,  the 
effects  of  the  pressure  gradient  on  the  turbulent  boundary  layer  structure  could  be  distinguished 
from  the  flat  plate  effects  at  the  same  location.  Figure  1-1  presents  a  schlieren  photograph  of  the 
present  adverse  pressure  gradient  flow  field.  This  figure  is  included  to  describe  the  present 
flowfield,  and  the  detailed  results  of  Figure  1-1  will  be  discussed  in  Chapter  7.  Clearly,  the 
shock  wave  has  coalesced  outside  of  the  boundary  layer.  The  approximate  turbulent  boundary 
layer  edge  has  been  designated.  The  weak  lip  shocks  are  results  of  misalignments  between  the 
test  sections. 

The  first  step  in  the  analysis  was  an  investigation  of  the  terms  that  require  modeling. 
Chapter  2  presents  the  governing  equations  of  turbulent  flow  and  a  discussion  of  the  utility  of 
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Flow  Direction 


Coalesced  Shock  Wave 


Weak  Step  Shocks 


Boundary  Layer  Edge 


Coalescing  Compression  Waves 


Figure  1-1 :  Flow  Over  Compression  Ramp  Model 


LDV  as  a  flowfield  analysis  tool.  Chapter  3  summarizes  the  key  principles  of  Laser  Doppler 


Velocimetry.  Thirdly,  the  experiment  was  designed  to  meet  the  objectives  described  in  Section 


1.3.  Chapters  4  and  5  present  the  test  facilities  and  equipment,  and  the  experimental  procedures 


developed  for  use  in  this  study,  respectively.  Chapter  6  presents  the  data  reduction  techniques 


used  to  extract  the  mean  and  turbulent  flow  properties.  After  the  research  scope  was  defined,  the 


testing  was  performed,  and  the  results  analyzed.  The  results  and  discussion  are  given  in  Chapter 


7.  Finally,  conclusions  and  recommendations  are  given  in  Chapter  8. 


Measurement  Locations 
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2.  Turbulence  Analysis 


Turbulent  flow,  which  is  governed  by  the  unsteady  Navier-Stokes  equations,  is 
characterized  by  temporal  and  spatial  scales  which  vary  over  several  orders  of  magnitude.  Thus, 
closed  form  solutions  to  practical  high  Reynolds  number  problems  are  not  possible  with  current 
or  near  future  computer  capabilities.  Hence,  engineers  and  scientists  must  rely  on  approximate 
time-average  methods  to  obtain  solutions  for  turbulent  flow.  Due  to  the  nonlinear  nature  of  the 
Navier-Stokes  equations,  this  procedure  results  in  additional  fluctuation  cross-correlation  terms. 
The  turbulence  “closure”  problem  is  the  science  or  art  of  developing  formulations  to  close  the 
system  (i.e.,  provide  enough  equations  for  the  unknowns).  This  chapter  presents  the  formal 
derivations  of  both  the  Reynolds  and  Favre  averaged  compressible  Navier-Stokes  equations. 


2. 1  Governing  Equations  ( Navier-Stokes  Equations ) 

The  compressible  Navier-Stokes  equations  govern  high  speed  flow  fields  where 
compressibility  effects  must  be  considered  [34],  Neglecting  body  forces,  the  continuity, 
momentum,  and  energy  equations  are  given,  respectively,  by 


3P  ,  ^(P“>)  _0 

dt  dxj 


(2.1) 


acpup  ^  3(pM,«y)  _  dp 

dt  dxj  dXj  dxj 
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where  the  stagnation  enthalpy,  shear  stress,  and  heat  flux  are  respectively  defined  as: 


h„  =h  +  —  U;U. 
"  2 


(2.4) 
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2.2  Reynolds-Averaged  Navier-Stokes  (RANS)  Equations 

The  Reynolds-averaged  Navier-Stokes  equations  [  1  ]  are  derived  when  time  averaged 
flow  properties,  like  the  velocity,  are  decomposed  as 


u{x,y,z,t)  =  u{x,y,z))+u'{x,y,z,t)  (2.7) 

where  m  is  the  time  averaged  velocity,  and  u  is  the  instantaneous  fluctuating  turbulent 
component.  The  time  averaged  quantity  M  is  defined  by 


M 


1 


u  dt 


(2.8) 


where  the  Af  must  be  finite  and  large  compared  to  the  period  of  random  turbulent  fluctuations, 
but  small  compared  to  the  mean  flow  time  scale. 

Once  the  decoupled  variables  are  substituted  into  the  Navier-Stokes  equations,  and  the 
rules  of  time,  averaging  applied  [1];  the  RANS  equations  for  continuity,  momentum,  and  energy 
are  obtained: 
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(2.9) 


ap  ^  3(pu^  +m;)_^ 
dt  dxj 

d  (pM,  +  p  x)  ^  3  {pUjUj  )  _  ap  ^  5  (^,7  +  X  .y  ) 
dt  dXj  dx-  dxj 

3(pg„  +PX)  ^  _  ^(ujXy+u'x'j -q. -q^) 

dt  dXj  dXj 

where  the  Reynolds  averaged  compressible  turbulence  terms  are  given  as 


m"  =  -p  u' 


t"  = 
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The  averaging  process  eliminates  several  key  turbulent  characteristics,  e.g.  frequency, 
phase,  and  wavelength  of  turbulent  motion  [3].  For  incompressible  flows,  the  apparent  mass  flux 
{m^ )  and  the  last  three  terms  of  the  turbulent  shear  stress  (Ty)  and  turbulent  heat  flux  (q/^ )  are 


zero. 


2-3 


2.3  Favre-Averaged  Navier-Stokes  (FANS)  Equations 

The  Favre-averaged  Navier-Stokes  equations  [3]  are  based  on  mass-weighted  time 
averages.  The  mass-weighted  time  average  velocity  is  decomposed  are  as 

u{x,y,z,t)  =  u(x,y,z)+u"{x,y,z,t)  ( 2.16) 

where  the  Favre  mean  value  is  defined  as: 


u 


pu 

P 


(2.17) 


and  u"  is  the  instantaneous  fluctuating  turbulent  component.  Note,  unlike  u  ,  m"  0.  The 
Favre-averaged  continuity,  momentum,  and  energy  equations  are: 


3p  ,  ^(p^y)_o 

dt  dXj 


(2.18) 
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where  the  fluctuating  Favre  averaged  turbulent  terms  are 


X^=-pu'UJ  (2.21) 


qf  =  phX 


and,  neglecting  viscosity  fluctuations. 


( 2.22) 
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( 2.23) 
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2.4  RANS  &  FANS  Relationship 

Because  of  the  nature  of  the  Favre-averaged  method,  there  is  no  explicit  appearance  of 
the  density-velocity  correlations  as  with  the  Reynolds-averaged  equations.  Thus,  the 
compressible  FANS  equations  closely  resemble  the  incompressible  Reynolds-averaged  form. 
For  this  reason,  Favre  averaging  is  almost  universally  adopted  in  the  aerospace  community  [3]. 

By  decomposing  a  flow  parameter,  say  the  streamwise  velocity  component,  u,  by  both 
methods,  and  equating,  the  relation  between  Reynolds  and  Favre-averaged  variables  is 
determined: 


Also, 


u  —  u 


=  u 


( 2.24) 


u"  =  u'  +  u"  (2.25) 

In  addition,  through  the  above  identities,  the  Favre  and  Reynolds  shear  stresses  are  related  by: 

pu'ii''  =  pu'u'j  +  pu"u'j  (  2.26) 

Since  the  second  term  on  the  right  hand  side  of  Eq.  (2.26)  is  fourth  order,  it  is  small  compared  to 
the  second  order  shear  stresses.  Thus, 


Therefore,  LDV  can  provide  accurate  estimates  of  the  Favre  averaged  shear  stress. 


( 2.27) 
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3.  Principles  of  Laser  Doppler  Velocimetry 


This  chapter  presents  a  brief  overview  of  the  basic  principles  behind  Laser  Doppler 
Velocimetry.  This  information  provides  a  sound  foundation  upon  which  the  proper  system 
settings,  given  in  Chapter  5,  were  based. 

3. 1  Overview 

The  Laser  Doppler  Velocimetry  (LDV)  system  has  been  developed  into  a  research  tool 
of  unique  capabilities.  This  technique  can  measure  particle  velocities  in  a  small  flow  volume 
that  are  very  close  to  a  solid  surface,  without  disturbing  the  flow  field.  As  the  name  implies, 
LDV  makes  use  of  lasers  and  a  Doppler  frequency  shift  to  calculate  the  velocity  of  a  gaseous 
flow.  Seed  particles,  such  as  smoke,  latex  spheres,  or  olive  oil,  are  introduced  into  the  flow  field. 
As  these  particles  pass  through  the  probe  volume,  the  laser  light  is  scattered.  Photomultiplier 
optics  are  used  to  measure  the  Doppler  shift,  which  is  the  difference  between  the  frequency  of 
the  scattered  and  transmitted  light.  Based  on  the  wavelength  and  optical  setup  of  the  laser,  the 
particle’s  velocity  can  be  measured  [10,  17]. 

A  3D  component  LDV  system  enables  velocity  measurements  in  three  orthogonal 
directions.  The  present  study  used  the  system  in  a  configuration  to  measure  velocity  components 
in  two  of  these  directions.  Measurements  in  each  direction  are  facilitated  by  a  single  channel, 
fringe  mode  setup.  The  combined  efforts  of  two  of  these  components  allowed  measurements  of 
the  u-  and  v- velocities  along  the  x-  and  y-directions,  respectively  (see  Figure  3-1). 
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3.2  Probe  Volume  &  Interference  Pattern 

The  transmitting  optics  split  each  wavelength  of  laser  light  into  two  beams.  These 
beams  are  passed  through  a  lens  which  causes  them  to  intersect.  The  intersection  of  the  beams 
creates  the  probe  volume  or  measurement  volume,  which  is  ellipsoidal  in  shape.  This  shape  is  a 
result  of  the  Gaussian  intensity  profiles  of  the  laser  beams  [8].  The  size  of  the  probe  volume  is 
determined  by  the  wavelength  and  diameter  of  the  laser  beam,  and  separation  distance  of  the 
emitted  beams.  Figure  3-1  shows  the  probe  volume. 


Figure  3-1:  Fringe  Disks  of  Probe  Volume  [8] 

The  intersection  of  the  two  coherent  laser  beams  creates  an  interference  pattern,  known 
as  interference  fringes.  Due  to  the  3D  nature  of  the  laser  beam,  the  interference  fringes  actually 
appear  as  fringe  disks.  Because  the  beams  intersect  at  the  beam  waists,  or  narrowest  width  of  the 
beam,  the  wave  fronts  are  nearly  parallel;  thus,  the  interference  fringes  will  be  parallel.  Parallel 
fringes  are  important  characteristics  because  of  the  resulting  uniform  velocity-frequency 
correlation  over  the  entire  probe  volume  [8].  The  fringe  spacing,  dj  ,  is  given  by 

,  V2 
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where  A,  is  the  laser  wavelength,  and  0  is  the  angle  between  the  two  incident  beams.  With  the 
fringe  spacing  determined,  the  velocity  of  the  particle  can  be  determined  by  measuring  the  time 
that  elapses  as  the  particle  passes  through  the  probe  volume.  Figure  3-2  shows  the  geometry  of 
the  intersecting  laser  beams. 


Direction  of  velocity  component,  u,. 


3.3  Calculation  of  Velocity 

The  laser  optics  include  both  transmitting  and  receiving  optics.  A  photomultiplier,  or 
PM  tube,  is  used  to  help  the  photodetector  sensor  convert  the  scattered  light  into  a  continuous 
signal.  The  frequency  of  the  signal  depends  on  the  motion  of  the  particle;  higher  frequencies  are 
associated  with  faster  moving  particles.  As  the  particle  traverses  the  probe  volume  and  crosses 
the  fringe  patterns,  a  series  of  “on”  and  “off’  signals  are  created.  From  this  signal,  the  time  to 
traverse  the  probe  volume  can  be  measured.  Since  both  the  fringe  spacing  and  traverse  time  are 
known,  the  particle  velocity  can  be  determined. 

The  direction  that  the  velocity  is  to  be  measured  is  assumed  perpendicular  to  the  bisector 
of  the  beam  intersection  angle  9  and  lies  in  the  plane  formed  by  the  transmitted  laser  beams.  The 
velocity  of  a  particle  passing  through  the  probe  volume  is  given  by 
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sin(0/2) 


(3.2) 


where  is  the  Doppler  frequency.  The  Doppler  frequency  is  determined  by  analyzing  the 
signal,  or  Doppler  burst,  created  as  a  particle  passes  through  the  probe  volume.  Figure  3-3  shows 
a  typical  Doppler  burst.  The  maximum  points  indicate  when  the  particle  has  passed  through  a 
fringe,  while  the  minimum  points  indicate  when  the  particle  is  between  fringes.  The  time 
interval  between  each  peak,  or  period  7^ ,  is  the  time  required  for  a  particle  to  cross  a  fringe  disk 
in  the  probe  volume.  The  corresponding  Doppler  frequency  is  given  by  =  1/7^  . 


Figure  3-3:  Doppler  Burst  [7] 

The  frequency  of  the  Doppler  shift  signal  is  determined  by  the  Burst  Spectrum  Analyzer 
(BSA).  The  BSA  breaks  the  signal  into  the  envelope  and  pedestal,  and  creates  a  Fast-Fourier 
Transform  (FFT)  to  match  the  Doppler  signal.  This  begins  when  the  signal  exceeds  a  threshold 
voltage,  and  continues  until  the  signal  drops  below  the  threshold  voltage.  The  length  of  time  that 
the  threshold  has  been  exceeded,  ,  is  recorded.  The  Doppler  burst  is  modeled  by  using  the 
record  length,  or  N  samples  taken  seconds  apart.  The  resulting  Doppler  frequency  is  given 
by 


(3.3) 
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If  the  envelope  of  the  Doppler  signal  exceeds  a  higher  threshold,  the  input  sampling  is  restarted. 
The  samples  stored  in  the  input  buffer  are  deleted,  and  the  sampling  process  is  repeated.  This 
may  be  repeated  several  times;  each  time  with  a  higher  threshold  voltage.  This  procedure 
ensures  that  the  Doppler  frequency  is  determined  by  the  part  of  the  burst  with  the  largest 
amplitude.  With  the  Doppler  frequency  determined,  the  particle’s  velocity  is  easily  determined 
by  Eq.  (3.2). 

The  bandwidth,  BW  ,  is  defined  as  the  range  of  frequencies  (or  velocities,  depending  on 
the  desired  units)  over  which  possible  frequencies  are  to  be  determined.  The  size  of  the 
bandwidth  limits  the  sample  time,  ,  and  its  associated  frequency,  =  1/T^  .  The  frequency  is 

proportional  to  the  bandwidth,  and  given  by  =  n(^BW)  where  n  is  greater  than  one,  and  pre¬ 
determined  by  the  BSA  configuration.  Therefore,  the  total  sampling  time  is  increased  as  the 
record  length  is  increased.  The  sampling  time  must  be  less  than  the  time  to  traverse  the  control 
volume  [7]. 

Each  of  two  components  of  velocity  are  measured  with  two  separate  BSA’s  and  laser 
beams  of  different  wavelength.  The  wavelength  of  each  laser  beam  is  programmed  into  the  BSA 
so  that  it  is  dedicated  to  measuring  velocities  in  the  plane  of  the  corresponding  laser  beams. 

Thus,  a  BSA  is  required  for  each  component  of  velocity  to  be  measured.  A  calibration  factor, 

,  gives  the  relationship  between  the  signal  frequency  and  velocity.  The  calibration  factor  is  a 

function  of  the  wavelength  and  optical  configuration,  and  is  given  as 


.  X/2 

^  sin(9/2) 


(3.4) 
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where  the  units  of  the  calibration  factor  are  /n/5/  MHz  .  The  calibration  factor  also  gives  the 
fringe  spacing,  ,  in  units  of  |J.m. 

3.4  Overcoming  Directional  Ambiguity 

The  problem  of  directional  ambiguity  exists  because  the  frequency  of  the  Doppler  signal 
is  independent  of  the  sign  of  the  direction  of  velocity.  The  same  Doppler  signal  would  be 
generated  for  positive  and  negative  velocities  of  the  same  magnitude.  A  phase  shift  is  introduced 
by  passing  one  of  the  beams  of  each  wavelength  through  the  Bragg  cell.  The  Bragg  cell  is  a  glass 
cell  that  includes  an  oscillating  electromechanical  transducer.  As  the  each  pair  of  beams  passes 
through  the  Bragg  cell,  a  frequency  phase  shift  is  introduced.  The  shift  in  frequency,  , 
causes  the  fringes  to  sweep  across  the  probe  volume  at  a  constant  velocity  in  the  direction  of  the 
desired  velocity  component.  Johnson  [16]  gives  this  constant  velocity  as 

Vf=d^f,  (3.5) 

where  is  the  fringe  spacing.  The  direction  that  the  particle  travels,  and  thus  its  corresponding 
velocity,  can  be  determined  from  the  sweeping  fringe  motion. 

3.5  Validation  of  Particles 

3.5.1  Unfortunately,  not  all  particles  that  traverse  the  probe  volume  give  accurate 
velocity  measurements.  Some  particles  pass  through  the  middle  of  the  volume,  while  others  just 
barely  clip  the  edge.  There  are  several  methods  employed  by  the  BSA  to  filter  “bad”  particles 
crossing  the  probe  volume.  Those  particles  that  pass  the  filtering  criteria  are  considered  valid 
particles. 
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3.5.1  Fringe  Count 

A  minimum  number  of  fringes  must  be  crossed  by  each  particle  before  it  can  be 
validated.  When  measuring  highly  turbulent  flows,  particles  may  enter  the  probe  volume  from 
any  direction.  Those  particles  entering  the  control  volume  nearly  parallel  to  the  fringes  will 
cross  fewer  fringes  than  particles  entering  perpendicular  to  the  fringes.  Because  the  desired 
measurement  direction  of  velocity  is  perpendicular  to  the  fringes,  particles  crossing  a  few  fringes 
nearly  parallel  are  not  accurately  representing  the  correct  velocity  component.  Also,  particles 
passing  perpendicular  to  the  fringes,  but  near  the  edge  of  the  probe  volume,  will  not  pass  through 
enough  fringes. 

3.5.2  Particle  Acceleration 

Particles  that  change  velocity  while  passing  through  the  probe  volume  are  not  validated. 
The  particle’s  acceleration  may  be  detected  by  comparing  the  time  to  cross  different  number  of 
fringes.  If  the  time  between  crossing  each  of  10  fringes  is  different  than  that  of  20  fringes,  then 
the  particle’s  velocity  is  not  constant.  Therefore,  the  accelerating  or  decelerating  particle  is  not 
valid,  and  its  contribution  is  rejected. 

3.5.3  Spectrum  Validation 

While  computing  the  Doppler  frequency,  the  BSA  compares  the  intensity  of  the  gathered 
light  at  each  sample.  For  example.  Figure  3-4  shows  the  intensity  at  32  different  sample 
locations.  The  two  largest  local  maximum  intensities  are  compared  to  the  overall  maximum 
intensity.  Particles  are  validated  if  the  local  maximum  is  four  times  less  than  the  global 
maximum  of  the  spectrum  [7]. 
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Figure  3-4:  Spectrum  Validation  [7] 

3.5.4  Doppler  Signal  Validation 

Another  form  of  validation  involves  analyzing  and  comparing  the  pedestal  and  envelope 
of  the  Doppler  burst.  The  pedestal  and  envelope  of  the  signal  are  generated  through  signal 
processing,  as  shown  in  Figure  3-5.  These  characteristics  of  the  Doppler  signal  are  compared  for 
particle  validation.  The  quality  factor,  which  is  a  user  defined  parameter,  is  the  ratio  between  the 
envelope  and  pedestal  peak  voltages.  Particles  must  have  signal  amplitudes  that  exceed  the 
quality  factor  to  be  validated. 
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Figure  3-5:  Envelope  and  Pedestal  Generation  [7] 

3.6  Coincidence  Window 

With  a  two-component  LDV  system,  a  particle  traversing  the  probe  volume  will  produce 
two  different  Doppler  signals;  one  corresponding  to  each  wavelength  of  light.  Because  the  two 
signals  are  independent  of  one  another,  there  must  be  a  method  of  correlating  the  two  bursts. 
Coincidence  processing  is  the  procedure  of  comparing  the  arrival  time  of  a  particle.  It  is  very 
important  that  the  BSA’s  are  tracking  the  same  particle  through  the  probe  volume  space. 
Otherwise,  the  two  signals  would  not  be  correlated,  the  velocity  correlations  and  turbulent  shear 
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stresses  would  be  artificially  low.  During  data  acquisition,  the  velocity,  arrival  time,  and  transit 
time  for  each  burst  are  recorded  by  each  BSA.  As  a  particle  enters  the  probe  volume,  the  first 
BSA  signals  the  second  BSA  to  “look”  for  the  same  particle.  The  coincidence  window  is  the 
maximum  allowable  difference  in  arrival  times  for  the  two  BSA  signals.  Coincident  bursts  are 
those  bursts  that  meet  this  criteria.  If  the  coincidence  window  is  exceeded,  then  the  BSA’s  are 
not  tracking  the  same  particle  and  the  burst  is  rejected.  The  size  of  the  coincidence  window  must 
be  carefully  selected  so  that  two  burst  created  by  a  single  particle  are  correlated,  yet  must  not 
allow  more  than  one  burst  from  one  BSA  correlate  with  a  burst  from  another  BSA.  Also,  the 
coincidence  window  must  be  a  fraction  of  the  time  that  a  particle  takes  to  traverse  the  probe 
volume. 

3. 7  Forward-  vs.  Back-Scattering  Method 

Doppler  bursts  are  created  by  collecting  and  analyzing  the  light  that  scatters  forward,  or 
ahead,  of  the  particle.  It  is  very  important  to  position  the  receiving  optics  in  a  configuration  so 
that  the  maximum  intensity  of  the  scattered  light  may  be  collected.  The  relative  magnitude  of  the 
scattered  light  depends  on  the  size  of  the  seed  particle,  as  seen  in  Figure  3-6.  When  the  particle 
diameter,  ,  is  larger  than  the  wavelength  of  the  laser  light,  a  larger  intensity  of  the  light  is 

scattered  ahead  of  the  particle.  Figure  3-7,  which  is  on  a  radial  log  scale,  indicates  that  the 
magnitude  of  the  forward  scattered  light  is  several  orders  of  magnitude  larger  that  backward 
scattered  light. 


dp<X. 

Figure  3-6;  Scattered  Light  Intensity  [6] 
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Figure  3-7;  Magnitude  of  Light  Intensity  [10] 

The  LDV  system  may  be  configured  in  either  back-  or  forward-scattering  collection 
mode.  In  order  to  take  advantage  of  the  forward  scattered  light,  the  receiving  and  transmitting 
optics  are  configured  in  forward-scattering  mode.  In  this  arrangement,  the  receiving  and 
transmitting  optics  are  on  opposite  sides  of  the  desired  measurement  region,  in  this  case  the  wind 
tunnel.  Forward-scattering  mode  enables  collection  of  the  larger  intensities  of  light,  resulting  in 
stronger  Doppler  signals. 

In  some  cases,  the  back-scattering  method  may  be  used.  In  this  arrangement,  the 
transmitting  optics  also  act  as  the  receiving  optics;  thus,  only  one  probe  is  required.  For  the 
present  system,  the  back-scattering  mode  is  a  “point-and-shoot”  method  because  the  optics  are 
pre-aligned  within  the  probe.  Accurate  compressible  measurements  are  typically  difficult  to 
obtain  due  to  the  low  power  of  the  scattered  light.  For  this  reason,  the  back-scattering  mode  was 
not  used  in  this  experiment. 

3.8  Seeding  Considerations 

The  quality  of  the  LDV  measurements  depends  greatly  on  the  method  of  seeding  the 
flow.  While  at  times  it  may  be  possible  to  make  accurate  LDV  measurements  from  particles  that 
occur  naturally  in  the  flow  field,  it  is  desirable  to  fdter  out  naturally  occurring  particles,  and  then 
seed  the  flow  with  particles  of  known  size  and  material  properties  [11].  Particles  must  be  small 
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enough  that  they  accurately  follow  the  fluid  flow  field,  yet  large  enough  to  scatter  enough  laser 
light  to  produce  a  measurable  Doppler  signal. 

Seeding  of  particles  can  occur  at  different  locations  in  the  flow  field.  In  Ref.  [11],  Elena 
makes  several  suggestions  for  seeding  supersonic  flows.  The  flow  field  upstream  of  the 
measurement  location  must  remain  undisturbed  by  introduction  of  the  particles.  Injecting  the 
particles  at  the  local  mean  flow  velocity  achieves  this  requirement.  Elena  also  showed  that  the 
wake  created  by  injecting  seed  particles  at  the  wall  greatly  affects  the  accuracy  of  boundary  layer 
measurements. 

Particle  sizing  and  material  composition  are  other  important  factors  that  affect  the 
accuracy  of  LDV  measurements.  Particles  must  be  small  enough  to  follow  even  the  smallest 
fluid  motions.  Particles  which  are  too  large  will  lag  behind  the  actual  flow  velocity,  thus 
experiencing  velocity  slip.  They  must  be  made  of  material  that  does  not  adhere  to  the  wind 
tunnel  or  windows.  When  shock  waves  are  studied,  the  particles  must  have  fast  response  times 
in  order  to  map  the  rapid  changes  in  flow  conditions.  High  density  particles  may  even  change  the 
flow  conditions  of  the  flow  field  being  analyzed. 

3.9  Optical  Alignment 

When  making  LDV  measurements,  it  is  possible  to  align  the  laser  beams  in  different 
ways  relative  to  the  desired  velocity  directions,  as  shown  in  Figure  3-8.  The  coordinate  system 
shows  the  direction  of  the  actual  flow  velocities,  while  the  black  and  white  circles  symbolize  the 
different  laser  beams.  The  lighter  arrows  indicate  the  direction  of  positive  velocity.  With  the 
first  method,  one  pair  of  laser  beams  measures  the  u-component  of  velocity,  while  the  v- 
component  is  measured  with  the  other  beams.  With  the  second  method,  the  beams  are  rotated  by 
an  angle,  in  this  case  45°  clockwise,  so  that  each  laser  is  measuring  a  contribution  of  each 
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component  of  velocity.  An  optical  transformation  must  be  used  during  post-processing  to  resolve 
the  components  of  velocity  in  the  actual  coordinate  system. 


Figure  3-8:  Possible  Optical  Alignments 

As  described  by  Elena  [11],  the  problem  of  angular  bias  is  introduced  when  the  optics 
are  aligned  as  on  the  right  side  of  Figure  3-8.  When  the  probe  volume  is  rotated,  the  number  of 
fringes  crossed  by  a  particle  could  be  reduced  as  the  particles  may  enter  the  volume  almost 
parallel  to  the  fringes.  As  mentioned  previously,  these  particles  are  not  validated.  Rejecting 
these  signals,  results  in  errors  in  the  calculation  of  turbulence  intensities  and  the  correlation 

u'v' .  Elena  and  LaCharme  [12]  reported  that  the  measured  values  of  turbulence  intensities  and 

magnitude  of  u'v'  are  decreased  when  the  optics  are  rotated  from  the  original  coordinate 
system. 
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4.  Test  Facilities  and  Equipment 


This  chapter  discusses  the  equipment  and  facilities  used  in  this  investigation.  The  entire 
experimental  investigation  was  conducted  in  the  AFIT  Mach  3.0  ( Re/m  =  1.75x10’ )  supersonic 
wind  tunnel.  This  chapter  describes  the  wind  tunnel  configuration  and  Laser  Doppler 
Velocimetry  system  used  in  this  experiment,  including  the  temperature  and  pressure  measuring 
apparatus,  and  flow  visualization  equipment. 

4.1  Wind  Tunnel  Facilities 

The  AFIT  supersonic  wind  tunnel  was  a  blow-down,  pressure- vacuum  system.  High- 
pressure  air  was  introduced  into  the  settling,  or  plenum  chamber.  The  end  of  the  tunnel  was 
connected  to  pressure- vacuum  tanks.  The  combination  of  high-pressure  air  and  a  vacuum 
provided  the  necessary  conditions  to  run  the  tunnel.  The  layout  of  the  AFIT  Mach  3.0  wind 
tunnel  is  presented  in  Figure  4-1. 


Flow  Straigtener  Nozzle  Gate  Valve 


Figure  4-1;  AFIT  Mach  3.0  Tunnel  Schematic 


4.1.1  Plenum  Chamber 


The  settling  chamber  pressure  and  temperature  were  maintained  at  2.06  ±  0.02  atm,  and 
295  ±  2  K  for  all  test  runs.  Plenum  chamber  pressure  and  temperature  were  measured  between 
the  flow  straightener  and  the  nozzle.  The  Pitot  pressure  was  measured  with  an  Endevco  0.69 
MPa  pressure  transducer,  and  the  temperature  was  measured  with  an  Omega  Engineering  Type  K 
thermocouple. 

4.1.2  Air  Supply 

Two  Atlas  Copco  GAU  807  air  compressors  provided  air  at  0.69  MPa.  The  air  was  dried 
by  two  Pioneer  Air  Systems,  Inc.  Model  R500A  Refrigerant  Air  Dryers  before  entering  the  high- 
pressure  reservoir.  The  mass  flow  rate  was  0.5  kg/s.  A  centrifugal  moisture  and  particle 
separator,  with  multi-layered  high  grade  filters,  was  used  to  cleanse  the  air  before  it  entered  the 
tunnel. 

Two  7.5  hp,  460  V  Stokes-Penwalt  Model  212-11  MicroVac  evacuation  pumps  were 
used  to  evacuate  the  16  cubic  meter  vacuum  tank  system.  The  vacuum  was  maintained  by  a 
Stokes-Penwalt  seal  exhaust  pump.  Each  pump  was  powered  by  a  Reliance  XW  Duty  Master 
AC  motor.  The  vacuum-pressure  system  empties  the  vacuum  tanks  to  2  mmHg,  or  267  Pa,  in 
approximately  7  minutes.  Wind  tunnel  test  runs  were  performed  when  the  pressure  in  the 
vacuum  tanks  was  below  15  mmHg,  or  2000  Pa.  This  configuration  provided  approximately  25 
seconds  of  constant  wind  tunnel  operation  before  unstarting  conditions  occurred. 

4.1.3  Supersonic  Nozzle 

The  converging-diverging  half-nozzle  generated  a  mean  flow  Mach  number  of  3.0.  The 
freestream  Reynolds  number  per  meter,  based  on  plenum  total  temperature  and  pressure,  was 
calculated  as  Re/m  =  1.75  x  10^ .  The  nozzle  exit  cross  section  was  6.35  x  6.35  cm,  and  the 
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distance  from  the  throat  to  the  end  of  the  nozzle  was  27  cm.  Figure  4-2  shows  the  Mach  3.0 
nozzle. 


Figure  4-2:  Mach  3.0  Nozzle 


4.1.4  Test  Section 

Two  test  sections  measuring  6.35  cm  by  6.35  cm  in  cross-section,  and  33.0  cm  in  length 
were  located  directly  downstream  of  the  nozzle.  The  test  sections  of  the  tunnel  were  designed  to 
facilitate  simple  interchange  of  models  for  other  users.  The  walls  of  the  test  section,  nozzle,  and 
diffuser  were  made  of  1.905  cm  thick  aluminum  alloy.  Except  for  the  flat  plate  measurements  at 
the  X  =  44  cm  location,  all  other  measurements  in  the  present  study  were  made  in  the  second  test 
section.  The  first  test  section  was  fitted  with  flat  plate  walls  on  all  sides. 

A  compression  ramp  model  was  fabricated  to  produce  an  adverse  pressure  gradient 
(APG).  After  evaluation  of  the  results,  it  was  determined  that  the  model  generated  regions  of 
both  adverse  and  favorable  pressure  gradient  (FPG).  The  FPG  was  generated  by  the  flow 
expanding  around  the  top  of  the  model.  The  model  was  carefully  designed  to  avoid  separated 
flows  in  the  regions  were  measurements  were  to  be  taken.  It  was  also  designed  so  that 
compression  waves  would  coalesce  into  a  shock  outside  of  the  boundary  layer.  The  curved  wall 
generated  a  mild  compression,  where  the  Clauser  pressure  gradient  parameter,  p  ,  was  roughly 
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1.12  at  68  cm  in  the  adverse  pressure  gradient  region,  and  approximately  -0.94  in  the  favorable 
pressure  gradient  region  at  7 1  cm.  The  compression  ramp  model  matched  the  curve  fit  given  by 

y(;c)=1.1858-0541x-l-0.0748x:"-0.0028x'  (4.1) 

where  F(x)  is  the  change  in  model  height  from  flat  plate  entrance  conditions,  and  the  x- 

coordinate  is  measured  from  the  front  of  the  test  section.  The  units  of  Eq.  (4.1)  are  centimeters. 
The  standard  deviation  of  the  curve  fit  was  0.0008  cm.  The  model  is  symmetric  about  the  mid¬ 
length.  Figure  4-3  shows  the  APG  model  along  with  the  possible  window  configurations  of  the 
test  section.  The  two  lines  represent  the  profiles  along  which  the  LDV  measurements  were 
made.  The  two  asterisks  on  the  profiles  indicate  the  location  of  the  origin  for  each  LDV 
measurement  profile.  The  model  was  mounted  as  the  ceiling  of  the  second  test  section. 


Figure  4-3:  Compression  Ramp  Model  with  Test  Section 

The  compression  ramp  model  also  included  pressure  taps  along  its  centerline.  The 
surface  pressure  was  measured  with  an  Endevco  0.103  MPa  pressure  transducer.  The  first  port 
was  located  3  cm  behind  the  front  of  the  model,  with  seven  more  ports  successively  spaced  2.5 
cm  apart.  A  surface  pressure  plot  is  included  in  Chapter  7. 

Because  a  forward-scattering  configuration  was  required,  each  side  wall  of  the  test 
section  was  modified  by  adding  two  7.62  cm  diameter  windows.  The  circular  optical  quality 
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glass  windows  were  1.27  cm  thick.  They  were  secured  in  the  walls  with  pressure  gaskets  and 
interchangeable  mounting  clamps.  The  moveable  clamps  allowed  adjustment  while  viewing 
different  sections  of  the  flow  field.  The  windows  were  offset  from  the  lateral  and  longitudinal 
centerline  of  the  wall.  Most  of  the  test  section  was  unobstructed  because  the  walls  could  be 
mounted  on  either  side  of  the  tunnel  test  section.  Figure  4-4  shows  the  possible  glass  window 
configurations  which  provided  excellent  observation  for  LDV  measurements  and  flow 
visualization  photographs. 


Figure  4-4:  Possible  Window  Configurations 

4.2  Laser  Doppler  Velocimetry  System 

The  Dantec  FiberFlow  three-component  Laser  Doppler  Velocimetry  system  used  in  this 
experiment  consisted  of  the  laser  and  necessary  optics  to  manipulate  the  beams,  and  a  three-axis 
traverse  that  allowed  precision  positioning  of  the  probe  volume.  The  data  processing  system 
included  three  Burst  Spectrum  Analyzers  (BSA),  and  software  to  facilitate  data  processing  and 
communication  with  a  PC  computer.  A  flow  seeder,  including  injector,  was  added  to  the  tunnel 
to  introduce  seed  particles  for  the  LDV  measurements. 

4.2.1  FiberFlow  System 

The  Dantec  Measurement  Technology  60X  FiberFlow  system  was  used  to  make  LDA 
measurements.  The  purpose  of  the  FiberFlow  system  was  to  generate  and  manipulate  the  laser 
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beams  necessary  for  LDV  measurements,  and  then  detect  and  record  the  resulting  Doppler  burst 
signals.  The  FiberFlow  system  consisted  of  a  laser  light  source,  transmitter,  photo  multiplier 
tubes,  and  optical  probes  with  lenses. 

The  Ion  Laser  Technology,  Model  5500A-(X)  Ar-ion  laser  provided  a  coherent  light 
source.  The  300  mW  laser,  capable  of  generating  laser  radiation  in  the  457  nm  to  514.5  nm  range 
was  optimized  to  275  mW.  The  laser  power  supply.  Ion  Laser  Technology  Model  5405A 
operating  at  210  Volts  and  20  Amps  required  a  220  VAC  power  source.  The  laser  was  designed 
to  operate  in  the  TEMoo  transverse  mode,  where  the  laser  power  is  centered  around  the  optical 
axis  with  maximum  power  in  the  center.  The  beam  diameter  of  the  laser  was  0.82  mm 
[16]. 

The  FiberFlow  60X41  transmitter,  which  was  equipped  with  six  60X24  manipulators, 
served  the  role  of  splitting  the  laser  light  source  into  three  separate  wavelengths,  and  then 
delivering  these  laser  beams  to  the  55X  optical  probes  [6].  The  laser  was  separated  into 
wavelengths  of  514.4  nm,  488  nm,  and  476.5  nm,  which  corresponded  to  the  green,  blue,  and 
violet  beams,  respectively.  Upon  entering  the  transmitter,  the  laser  light  passed  through  the 
Bragg  cell  which  introduced  a  +40  MHz  frequency  shift,  fg  ,  to  one  beam  of  each  wavelength. 

After  passing  through  the  Bragg  cell,  the  beams  passed  through  the  six  60X24 
manipulator.  The  manipulators  positioned  the  laser  beam  on  the  exact  center  of  the  fiber  optic 
cable  by  causing  the  beams  to  converge.  Each  manipulator  has  four  adjustments  which 
controlled  the  beam  angle  and  displacement  position  of  the  laser  beam  relative  to  the  optical  axis 
of  the  fiber  optic  cable.  The  path  of  the  laser  beams  ran  from  the  manipulators  through  the  fiber 
optic  cables  to  the  optical  probes  [8].  Figure  4-5  shows  the  FiberFlow  system. 
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Figure  4-5:  Laser  Source  and  FiberFlow  System 
Each  of  the  60X  optical  probes  included  a  fiber  optic  transducer,  distributor  unit,  and 
transmitting  and  receiving  fibers  [8].  The  distributor  unit  served  as  an  interface  between  the 
transducer  head,  which  was  connected  to  the  fiber  cable,  and  the  transmitting  and  receiving 
fibers.  Two  transmitting  fibers  are  included  for  each  laser  beam  component;  the  2D  probe  had 
four  transmitting  fibers.  Two  pair  of  beams,  the  blue  and  green  beams,  ran  through  the  60X61 
2D  probe,  while  the  remaining  violet  beam  ran  through  the  ID  60X60  probe.  A  55X12  beam 
expander,  which  was  used  to  reduce  the  measurement  volume,  was  mounted  on  the  ID  probe. 

The  receiving  fibers  returned  the  acquired  signal  to  the  photo  multipliers.  The  focal  length  of  the 
60  mm  diameter  lens  on  each  probe  was  600  mm.  Each  probe  was  secured  to  a  micrometer 
adjustable  mount.  The  2D  probe  mount  was  only  able  to  swivel  in  the  horizontal  plane.  The  ID 
probe,  which  was  used  as  the  receiving  optics  in  the  forward-scattering  configuration,  was 
mounted  on  a  stand  that  included  horizontal  swivel  and  vertical  adjustments. 
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Figure  4-6:  ID  Probe  with  Adjustable  Base 


To  make  LDV  measurements,  the  lasers  beams  must  intersect  at  the  same  location. 

While  the  476.5  nm  wavelength  beam  did  not  actually  measure  a  velocity  component,  alignment 
with  the  probe  volume  of  the  2D  lasers  ensured  proper  reception  of  the  scattered  light.  Poor 
signal  quality  resulted  when  the  probe  volumes  of  the  transmitting  and  receiving  optics  were  not 
properly  aligned.  The  probe  mounts  were  positioned  on  the  ends  of  a  1.83  m  mounting  bench  on 
the  traverse.  The  angle  of  offset  for  each  probe  was  3.5°  from  the  mounting  bench  axis.  Based 
on  this  offset  angle  and  the  600  mm  focal  lengths  of  each  lens,  the  lenses  of  the  receiving  and 
transmitting  probes  were  separated  by  a  distance  of  1 .2  meters. 

Scattered  light  signals  entered  the  receiving  fibers  which  were  connected  to  the  photo 
multiplier  tubes  via  the  fiber  optic  cables.  When  high-voltage  was  applied  to  the  PM  tubes,  the 
scattered  light  signal  was  converted  into  a  Doppler  frequency  signal.  The  55X35  dual¬ 
component  tube  was  used  for  measurement  of  the  u-  and  v-components  of  velocity.  While  this 
Dantec  LDV  system  could  also  be  configured  for  three-component  LDV  measurements,  only  the 
dual-component  PM  tube  was  used  in  this  experiment. 


4.2.2  Dantec  3D  Traverse  System 

The  probe  volume  for  LDA  measurements  was  positioned  with  the  aid  of  the  Dantec  3D 
traverse.  Traverse  movement  for  each  axis  was  controlled  by  a  Whedco  DC  servo  motor  and 
encoder.  The  encoder  converted  1000  pulses  into  1  revolution  of  the  traverse  axis  load  screw. 
Each  revolution  of  the  load  screw  moved  the  traverse  2  mm.  Traverse  positioning  was  accurate 
to  80|im  over  a  600  mm  range.  The  three-axis  traverse  allowed  maximum  travel  of  600  mm  in 
each  of  three  axis  directions  [5].  Each  axis  was  fitted  with  limit  switches  to  prevent 
displacements  that  would  exceed  the  physical  limitations  of  the  traverse.  Figure  4-7  shows  the 
traverse  configuration  with  mounting  bench  and  optical  probes. 


Figure  4-7:  3D  Traverse  with  Probes 

The  y-axis  of  the  traverse  was  aligned  with  the  wind  tunnel,  with  the  z-axis  running 
perpendicular  to  the  horizon.  The  traverse’s  x-axis  completed  a  standard  right-hand-rule 
coordinate  system.  The  traverse  was  controlled  via  the  BURSTware  software.  The  mounting 
bench,  discussed  in  the  previous  section,  was  aligned  with  the  x-axis  of  the  traverse,  or 
perpendicular  to  the  wind  tunnel. 
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4.2.3  Burst  Spectrum  Analyzers 

The  purpose  of  the  BSA’s  was  to  perform  Fast  Fourier  transforms  (FFT)  on  the  LDA 
signal.  Three  BSA’s  were  included  in  the  Dantec  3D  Laser  Doppler  Velocimetry  system. 
Because  only  two  components  of  velocity  were  measured,  only  two  of  the  three  BSA’s  were 
used.  One  of  the  BSA’s  was  a  Dantec  57N20  BSA  Enhanced  model  which  included  keyboard 
and  LED  displays  on  its  front  panel.  The  remaining  two  analyzers  were  Dantec  57N35  BSA 
Enhanced  Model  S  units.  The  57N20  unit,  BSA  #1  which  contained  the  1  MHz  reference  master 
clock,  measured  the  velocity  component  which  corresponded  to  the  514.5  nm  wavelength  (green) 
beam.  The  second  BSA,  BSA  #2,  corresponded  to  the  488  nm  wavelength  (blue)  beam. 

Because  the  57N35  units  require  a  master  clock  for  coincidence  filtering,  they  were 
connected  to  the  57N20  unit  with  a  “sync.bus”  connector.  The  40  MHz  frequency  shift,  required 
by  the  Bragg  cell,  was  generated  by  BSA  #1.  The  cables  from  each  photo  multiplier  tube  were 
connected  to  the  BSA  corresponding  to  the  particular  wavelength  of  light.  The  BSA’s  were 
connected  to  the  PC  computer  via  an  IEEE  488  standard  connector  interface  [7]. 

4.2.4  PC  Computer  and  BURSTware  Software 

The  BSA  settings  were  controlled  with  Dantec  BURSTware  version  3.00  software, 
which  was  installed  on  the  local  PC  computer  [5].  This  program  served  as  an  interface  between 
the  user  and  3D  LDV  system.  The  computer  was  a  Gateway  2000  486DX/33  PC  with  4  Mbytes 
of  RAM  and  a  200  Mbyte  harddrive.  The  BURSTware  program  was  configured  to  interface  with 
all  three  BSA’s.  The  only  method  for  controlling  BSA  #2  and  #3,  because  they  lacked  display 
and  keyboard  controls,  was  with  this  software.  The  traverse  was  also  controlled  with  the 
software.  BSA  settings,  such  as  bandwidth,  center  frequency,  PM  tube  high-voltage,  and  data 
collection  mode,  could  all  be  controlled  and  varied  with  the  BURSTware  software.  The 
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processed  data  was  exported  from  the  software,  and  used  with  other  in-house  software  [23]  to 
present  turbulent  and  mean  flow  results. 


4. 3  Seeder  A pparatus 

The  seeder  apparatus  consisted  of  two  components;  the  atomizer  and  an  injector.  In 
order  to  introduce  particles  of  known  size  and  material  property  into  the  flow  field,  a  TSI 
Incorporated  Model  9306  6-jet  Atomizer  was  used.  The  seed  particles  passed  through  a  Tygon 
tubing  hose  placed  in  an  opening  in  the  settling  chamber.  An  injector  was  mounted 
approximately  4  cm  upstream  of  the  flow  straightener.  Figure  4-8  shows  the  atomizer  and  seeder 
configuration. 


Figure  4-8:  Six-jet  Atomizer  and  Injector  Schematic 

The  seed  particles  were  generated  by  passing  100%  pure  olive  oil  through  the  atomizer. 
Air  at  0.207  MPa  entered  the  atomizer  and  produced  an  aerosol  flow  rate  per  jet  of 
approximately  7450  cmYmin .  A  linear  relationship  between  flow  properties  and  the  number  of 
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jets  existed;  the  flow  rate  doubled  when  two  Jets  were  used.  The  diameter  of  the  olive  oil 
particles  was  nominally  0.6  |im[29]. 

The  flow  injector  was  made  of  1 .25  cm  diameter  copper  tubing.  A  90°  elbow  was  fitted 
on  one  end  of  the  injector.  The  copper  tubing  was  flattened  to  reduce  the  downstream  wake 
produced.  About  25  holes,  which  spanned  a  length  of  7.5  cm,  were  cut  in  the  top  edge  of  the 
injector,  which  was  mounted  vertically  in  front  of  the  flow  straightener..  The  Tygon  tubing 
passed  through  a  hole  in  the  settling  chamber  and  connected  to  the  atomizer  outlet. 

4.4  Flow  Visualization 

Schlieren  and  shadowgraph  techniques  were  employed  to  make  flow  visualization 
studies.  Figure  4-9  shows  the  flow  visualization  setup.  A  Polaroid  camera  using  Type  57  black 
and  white  Polaroid  film  was  used  along  with  mirrors  of  focal  length  1.53  m.  The  light  source 
was  a  Xenon  Nanopulser  Model  437B  with  N787B  spark  light.  The  flash  duration  was  10 
nanoseconds. 


Figure  4-9:  Shadowgraph  and  Schlieren  Setup  [15] 
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4.5  Computer  Resources 


4. 5. 1  Personal  Computer 

The  author’s  home  personal  computer  was  used  for  data  analysis,  and  document 
preparation.  The  computer  was  a  486DX33  PC  with  8  Mbytes  of  RAM  and  a  combined 
harddrive  space  of  550  Mbytes.  Software  resources  included  Microsoft  Office  v.4.3 
applications,  Mathsoft  MathCad  V.5.0+,  and  Golden  Software  Grapher  v.1.23.  Documents  were 
printed  on  a  Hewlett-Packard  DeskJet  500  printer.  Modem  connection  to  AFIT  provided  file 
transfer  and  remote  login  capability  to  the  Air  Force  Institute  of  Technology  UNIX  computer 
network. 

4.5.2  AFIT  Computer  Resources 

The  computer  resources  at  the  school  included  Sun  Sparcstation  20,  running  in  UNIX 
format,  and  laboratory  PC,  running  in  DOS  and  Windows.  Network  connections  provided 
support  for  file  transfer,  data  processing,  and  data  storage. 
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5.  Experimental  Procedures 


This  chapter  describes  the  experimental  procedures  used  for  the  present  experiment.  The 
primary  measurement  technique  was  Laser  Doppler  Velocimetry  (LDV).  As  mentioned  in 
Chapter  4,  the  LDV  system  contained  the  traverse,  optics,  and  software  necessary  to  perform  the 
experiments.  A  description  of  data  acquisition  is  also  included. 

5.1  Traverse  Positioning 

The  traverse,  described  in  Chapter  4,  was  placed  under  the  wind  tunnel,  which  was  about 
1.0  meter  above  the  floor,  with  its  y-axis  aligned  with  the  tunnel.  Thus,  moving  the  traverse 
along  its  y-axis  moved  the  measuring  volume  along  the  tunnel’s  centerline.  The  2D  optical  probe 
was  used  to  accurately  align  the  traverse.  The  probe  was  rotated  on  its  mounting  stand  until  the 
angular  variation  from  the  tixis  of  the  mounting  bench,  0  ,  was  zero.  Alignment  marks  on  the 
mounting  stand  aided  with  this  alignment.  With  the  probe  in  this  configuration,  the  laser  was 
turned  on  causing  the  beams  to  pass  through  the  optical  glass  windows  on  the  tunnel.  The  laser 
beams  produced  reflections  in  the  glass  windows.  The  traverse  was  shifted  until  the  reflections 
of  the  four  laser  beams  returned  into  the  2D  optical  probe.  The  axes  of  the  traverse  were 
considered  both  parallel  and  perpendicular  to  the  tunnel  when  this  alignment  was  achieved. 

5.2  Optical  Alignment 

Several  important  factors  were  considered  for  correct  optical  alignment.  These  factors 
included  forward-scatter  setup,  optical  axis  alignment,  measurement  profile  origin  positioning, 
and  laser  beam  intersection  adjustment. 
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5.2. 1  Forward-scatter  Mode 


The  Dantec  3D  Laser  Doppler  Velocimetry  system  can  be  configured  for  both  forward- 
and  back-scatter  collection  modes.  It  was  expected  and  verified  during  the  present  experiment, 
that  the  back-scatter  mode  would  not  generate  a  measurable  signal.  Hence,  the  receiving  and 
transmitting  optics  were  configured  for  forward-scatter  collection  mode,  as  directed  in  Ref.  [6] . 

In  forward-scattering  mode,  the  receiving  and  transmitting  optics  are  on  opposite  sides  of 
the  region  to  be  tinalyzed;  in  this  case,  the  test  section  of  the  wind  tunnel.  The  probe  mounts 
were  placed  near  each  end  of  the  traverse’s  mounting  bench.  The  offset  angle,  0  ,  between  the 
centerline  of  the  mounting  bench  and  the  direction  of  the  emitted  laser  beams  was  set  at  3.5°. 
Figure  5-1  shows  the  forward-scattering  setup.  Both  the  ID  and  2D  optical  probes  were  oriented 
at  this  angle.  This  was  the  minimum  angle  that  would  allow  the  lasers  from  the  transmitting 
probe  to  pass  through  the  tunnel  and  not  directly  intersect  the  lens  of  the  receiving  probe.  This 
setup  was  deemed  necessary  to  insure  that  the  very  sensitive  photo  multiplier  tubes  were  not 
destroyed.  The  bias  error  introduced  into  the  velocity  measurements  was  found  to  be  negligible 
(See  Appendix  A). 


Figure  5-1:  Forward-scattering  Configuration 
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5.2.2  Optical  Axis  Alignment 

As  mentioned  in  Chapter  3,  several  optical  axis  configurations  may  be  used  in  LDV 
measurements.  To  avoid  the  problems  associated  with  angular  bias  [1 1],  the  optical  axes  were 
aligned  to  directly  measure  the  u-  and  v-components  of  velocity.  This  process  also  negated  the 
requirement  of  providing  an  optical  transform  for  post-processing.  The  optical  axes  of  the  2D 
probe  were  aligned  with  model  test  section  body  intrinsic  coordinates;  the  u-component  of 
velocity  was  aligned  parallel  to  the  wall,  with  the  v-component  of  velocity  normal  to  the  wall.  In 
the  case  of  the  flat  plate,  the  u-component  of  velocity  was  aligned  with  the  tunnel’s  x-axis,  which 
corresponded  to  the  y-axis  of  the  traverse. 

The  orientation  of  the  optics  depended  on  the  curvature  of  the  model.  Using  Eq.  (4.1), 
the  slope  of  the  line  tangent  to  the  model  surface  at  a  particular  x-location  was  given  as 

y'(^)  =  -0-541  +  0.14956(x)  -  0.0084(x' )  (5.1) 

where  x  is  measured  in  centimeters  from  the  front  of  the  test  section.  LDV  measurements  were 
made  at  two  locations  for  the  adverse  pressure  gradient  case,  x  =  68  cm  and  x  =  71  cm  behind 
the  throat  of  the  nozzle,  corresponding  to  positions  8  and  1 1  cm  behind  the  front  of  the  model, 
respectively.  Flat  plate  measurements  were  also  made  at  x  =  71  cm,  for  comparison  with  the 
adverse  pressure  gradient  case.  These  locations  were  determined  by  raising  the  traverse  so  that 
the  probe  volume  was  positioned  just  above  the  tunnel.  An  index  card  was  placed  on  the 
centerline  of  the  tunnel,  and  the  location  of  the  beam  was  positioned  the  required  distance  behind 
the  front  of  the  model.  At  x  =  68  cm,  the  slope  of  the  tangent  line  was  0.1 1788.  At  x  =  71  cm, 
the  slope  of  the  tangent  line  was  0.08776.  These  slopes  of  the  tangent  lines  at  68  and  71  cm 
corresponded  to  angular  rotations  of  ([)  =  6.72  and  (|)  =  5.02  ,  respectively,  where  (\)  was 
measured  relative  to  the  horizon. 
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Once  the  proper  distance  from  the  nozzle  was  determined,  the  traverse  was  lowered  until 
the  beams  were  on  the  edge  of  the  model.  Clamps  on  the  probe  mounting  were  loosened,  and  the 
2D  optical  probe  was  rotated  through  the  angle  (()  until  the  488  nm  beam  (blue)  appeared  to  be 
aligned  parallel  to  the  wall.  Because  the  beams  emitted  from  the  probe  were  perpendicular,  the 
514.5  nm  beam  (green)  was  automatically  aligned  normal  to  the  wall.  In  this  configuration,  the 
blue  and  green  beams  were  set  to  measure  the  u-  and  v-components  of  velocity,  respectively. 
Slight  adjustments  in  the  height  of  the  traverse  were  required  during  this  alignment.  The 
orientation  of  the  beams  from  the  ID  receiving  optics  did  not  need  to  be  adjusted  relative  to  the 
2D  probe  beam  orientation. 

5.2.3  Establishment  of  Data  Profile  Origin 

Prior  to  making  LDV  measurements  across  the  boundary  layer,  the  origin  for  the  data 
profile  needed  to  be  established.  This  location  would  provide  a  reference  point  from  which  the 
probe  volume  was  adjusted.  With  the  probe  optics  already  aligned  with  the  model  coordinates, 
establishing  the  origin  was  achieved  by  adjusting  the  height  of  the  laser  beams  until  the  “4  spot 
pattern”  was  centered  on  the  edge  of  the  model.  In  particular,  the  traverse  was  positioned  so  that 
half  of  each  blue  beam  passed  through  the  tunnel,  and  the  other  half  was  clipped  by  the  model. 
Figure  5-2  shows  the  origin  configuration.  The  assumption  was  made  that  the  tunnel  was  leveled 
in  both  its  x-  and  y-axes. 


Data  Profile  Origin 


Figure  5-2:  Data  Profile  Origin 


5-4 


5.2.4  Laser  Intersection  Adjustment 

With  the  optical  axes  aligned  with  the  tunnel,  and  data  collection  origin  established,  the 
next  step  in  the  process  was  adjustment  of  the  laser  beams  to  ensure  that  the  receiving  optics 
were  focused  on  the  measurement  volume,  as  described  in  Chapter  3.  The  violet  (476.5  nm) 
components,  which  emitted  from  the  receiving  optics,  were  used  to  focus  the  optics.  The  2D 
transmitting  probe  was  adjusted  such  that  the  6  beams,  2  of  each  wavelength,  all  intersected  at 
the  same  region  in  space. 

Fine  tune  alignments  were  monitored  with  the  BURSTware  program.  There  were  two 
options  available  for  data  collection  mode;  online  display,  and  acquire  and  save.  These  modes 
determine  whether  the  acquired  data  is  saved  to  disk,  or  only  presented  on  the  monitor.  Online 
display  mode  was  used  to  adjust  the  lasers.  The  BURSTware  was  configured  so  that  each  BSA 
had  a  center  frequency  of  zero,  and  a  bandwidth  of  about  30  m/s.  (Details  of  the  BURSTware 
configuration  are  found  later  in  this  chapter.)  The  seeder  was  turned  on  so  that  olive  oil  particles 
were  introduced  into  the  tunnel.  Because  neither  the  high-pressure  air,  or  vacuum-suction  was 
applied  during  laser  alignment,  the  velocity  of  the  particles  was  nearly  zero.  With  the  laser 
turned  on,  the  particles  were  illuminated  by  the  laser.  The  scattered  light  Doppler  signals  were 
detected  by  each  BSA.  Data  were  acquired  for  approximately  5  seconds.  The  online  display 
showed  the  validation  percent  of  each  velocity  component  for  the  current  alignment.  By 
adjusting  the  micrometers  of  the  ID  probe,  the  beam  intersection  alignment  was  adjusted  until 
the  validation  percentages  for  both  components  were  above  90%.  The  validation  percent  was  the 
ratio  of  the  number  of  validated  bursts  to  the  total  number  of  detected  bursts.  The  data  rate  and 
number  of  bursts  collected  were  also  important  considerations  during  laser  alignment.  Data  rates 
and  number  of  collected  bursts  below  4  kHz  and  4000  bursts,  respectively,  were  unacceptable. 
Despite  a  high  validation  rate,  the  lasers  were  adjusted  if  either  of  these  parameters  was  low. 


5-5 


This  trial  and  error  task  proved  to  be  challenging  and  time  consuming.  The  validation 
percentages  were  very  sensitive  to  slight  adjustments;  in  fact,  different  validation  percentages 
resulted  when  the  micrometers  were  returned  to  previous  settings.  After  the  intersection  of  the 
laser  beams  was  optimized,  the  optics  of  the  LDV  system  were  ready  to  make  measurements. 

5.3  BURSTware  Configuration 

5.3.1  Parameter  Files 

The  profiles,  where  LDV  measurements  were  taken,  were  located  perpendicular  to  the 
wall  with  the  origin  at  the  previously  mentioned  locations.  Data  were  collected  at  about  15 
locations  across  the  boundary  layer.  The  BURSTware  software  required  a  parameter  file  for 
each  location.  The  parameter  file  included  all  the  traverse  and  BSA  settings  to  make  the  LDV 
measurements.  All  files  associated  with  the  same  profile  had  the  same  filename.  The  length  of 
these  root  names  was  limited  to  four  characters.  The  run  number,  which  corresponded  to  the 
location,  was  added  to  the  root  name.  Parameter  filenames  were  limited  to  six  characters 
because  the  software  used  the  last  two,  of  the  8  character  DOS  filename,  for  other  file 
designations  [5],  Examples  of  parameter  filenames  were  ADVlOl.par,  ADV102.par, 
FLATOl.par,  FLAT02.par,  etc.  The  parameter  files  at  each  location  were  identical,  except  for  a 
few  variations.  The  center  frequency,  which  corresponded  to  the  mean  velocity,  was  slightly 
increased  as  the  LDV  measurement  location  approached  the  edge  of  the  boundary  layer. 

5.3.2  BURSTware  Menu  Settings 

By  selecting  the  Setup  and  Acquire  option,  the  BURSTware  entered  the  mode  that  would 
allow  input  of  the  necessary  parameters.  This  option  enabled  parameter  file  management,  optical 
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setup,  BSA  and  traverse  control,  and  data  collection  parameter  assignment.  This  section  details 
the  BSA  settings  used  for  the  current  investigation. 

The  BURSTware  software  required  information  concerning  the  optics  used  by  the  LDV 
system.  These  settings  were  entered  with  the  Optics  option.  The  calibration  factor,  ,  which 

provided  a  direct  frequency- velocity  relationship,  was  based  on  the  laser  beam  wavelength,  beam 
separation  distance,  and  focal  length  of  the  lens  (see  Chapter  3).  The  calibration  factor  was 
equivalent  to  the  fringe  spacing  in  |lm.  Table  5- 1  summarizes  the  optical  settings  used  for  this 
investigation. 


Table  5-1:  BURSTware  Optical  Menu  Settings 


BSA#1 

BSA  #2 

Calibration 

(m/s/MHz) 

8.128 

7.709 

Wavelength 

(nm) 

514.5 

488.0 

Beam  Separation 

(mm) 

38 

38 

Focal  Length 

(mm) 

600 

600 

The  size  of  the  measurement  volume,  given  in  Table  5-2,  was  a  direct  result  of  these 
optical  settings.  The  coordinate  system  used  by  the  BURSTware  is  the  same  as  in  Figure  3-1. 
The  probe  volumes  for  both  sets  of  laser  beams  were  very  elongated. 


Table  5-2:  Control  Volume  Measurements  Settings 


BSA#1 

BSA  #2 

Gaussian  Beam  Diam. 

(mm) 

1.35 

1.35 

Beam  Expander  Ratio 

1.0 

1.0 

Beam  Collimator  Ratio 

1.0 

1.0 

AX 

(mm) 

0.291 

0.276 

AY 

(mm) 

0.291 

0.276 

AZ 

(mm) 

9.199 

8.725 
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When  the  Acquire  option  was  selected,  specific  data  acquisition  parameters  could  be 
entered  using  the  traverse  location,  Quick  setup.  Softkey  setup.  Program  setup,  and  Validation 
menu  options.  The  traverse  was  moved  to  the  proper  data  collection  position  before  other 
settings  were  adjusted.  This  precaution  was  taken  to  ensure  that  the  laser  beams  did  not  directly 
reflect  off  the  tunnel  and  back  into  the  optics.  Such  an  event  could  permanently  damage  the 
photo  multiplier  tubes.  Since  most  of  the  measurements  were  taken  in  the  center  of  the  tunnel, 
the  x-location  of  the  traverse  was  not  moved.  These  locations  originated  at  the  origin  and  were 
perpendicular  to  the  wall.  The  equations  for  the  profile  at  x  =  68  was  given  by 

Z  =  -8.4832y  (  5.2) 

and,  the  equation  for  the  profile  at  x  =  7 1  cm  was 

z  =  -11.3947y  (5.3) 

These  equations  were  derived  by  using  the  slopes  of  the  tangent  line  at  each  location  from 
Chapter  4.  The  slope  of  the  normal  line  was  the  negative  inverse  of  the  slope  of  the  tangent  line. 
A  few  traverse  coordinates  for  a  sample  profile  across  the  boundary  layer  at  x  =  68  cm  are  given 
in  Table  5-3.  While  the  accuracy  of  the  traverse  was  much  greater,  the  BURSTware  software 
limited  the  traverse  location  to  only  ±0.01  mm,  which  as  discussed  in  the  error  analysis  in 
Appendix  A,  was  the  largest  source  of  uncertainty  in  the  velocity  measurements  near  the  wall. 


Table  5-3:  Sample  Traverse  Coordinates 


Parameter  file  number 

Y  coordinate 

Z  coordinate 

01 

0.07 

-0.58 

02 

0.09 

-0.76 

03 

0.11 

-0.93 

04 

0.15 

-1.27 

05 

0.2 

-1.70 

06 

0.25 

-2.12 

07 

0.3 

-2.54 
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The  settings  of  both  BSA’s  were  controlled  from  the  Quick  menu.  These  settings  are 
detailed  in  Table  5-4.  The  center  frequency  and  bandwidth  gave  the  velocity  range  over  which 
possible  bursts  may  be  detected.  The  largest  possible  bandwidth  was  selected  to  ensure 
collection  of  a  wide  range  of  velocities.  With  the  frequency  shift  enabled,  a  wider  range  of 
center  frequency  and  bandwidth  settings  was  possible.  The  record  length  was  the  number  of 
samples  used  by  the  BSA  to  digitize  a  Doppler  burst.  Signal  gain  amplified  the  incoming  signal, 
while  the  pedestal  attenuation  reduced  the  pedestal  component  of  the  input  signal.  The 
maximum  anode  current  was  limited  by  the  type  of  photo  multiplier  tubes  [8].  The  high  voltage 
to  the  photo  multiplier  tubes  needed  to  be  enabled  to  allow  for  measurements  within  the  specified 
velocity  range.  The  high  voltage  to  BSA  #2  was  higher  because  the  power  of  the  green  beams 
was  less  than  that  of  the  blue  beams. 


Table  5-4:  BURSTware  Quick  Menu  Settings 


BSA#1 

BSA  #2 

Velocity  Units 

m/s 

m/s 

Center  Frequency 

u 

00 

o 

=  390 

Bandwidth 

325.12 

246.7 

Shifter  Mode 

Norm 

Norm 

Record  Length 

16 

16 

Signal  Gain 

(dB) 

42 

42 

High  Voltage 

(V) 

1304 

1504 

Max.  Anode  Current 

1.6 

1.6 

Pedestal  Attenuation 

6 

6 

Duty  Cycle 

(%) 

100 

100 

Dead  Time 

(msec) 

0 

0 

High  Voltage  On  ? 

Yes 

Yes 

After  the  Quick  menu  settings  were  entered,  the  Program  menu  settings  were  entered. 
The  number  of  bursts  to  be  collected,  and  type  of  data  to  be  recorded  were  selected  with  this 
option.  For  the  present  investigation,  data  was  collected  for  12  seconds,  or  until  150000  bursts 
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were  recorded  by  either  BSA.  Since  the  timeout  criteria  was  always  fulfilled  before  the 
maximum  number  of  bursts  was  collected,  the  number  of  bursts  varied  from  location  to  location. 
Near  the  wall,  the  number  of  bursts  was  lower  than  that  in  the  freestream.  In  any  event,  the 
minimum  number  of  bursts  at  a  single  location  was  6000.  The  remaining  three  settings 
determined  what  type  of  data  was  recorded  in  the  raw,  or  pre-coincidence  filtering,  data  files. 
Transit  time  is  the  amount  of  time  the  particle  was  in  the  control  volume.  Arrival  time,  which 
must  be  enabled  for  coincidence  filtering,  related  the  time  the  particle  was  first  detected  by  the 
BSA.  Obviously,  the  velocity  must  be  recorded  to  enable  analysis  of  velocity  fluctuations  and 
turbulent  quantities.  Table  5-5  shows  the  exact  settings  used  in  this  experiment. 


Table  5-5:  BURSTware  Program  Menu  Settings 


BSA#1 

BSA  #2 

Timeout  (sec) 

12 

12 

Number  of  Bursts 

150000 

Velocity 

Yes 

Yes 

Transit  Time 

Yes 

Yes 

Arrival  Time 

Yes 

Yes 

The  Softkey  menu  option  allowed  BSA  settings,  not  included  in  the  Quick  menu  option, 
to  be  entered.  Table  5-6  lists  the  settings  used  in  this  study.  BSA  #2  used  the  master  clock 
contained  in  BSA  #1.  The  coincidence  mode  settings  were  selected  so  that  the  second  BSA 
would  search  for  a  burst  after  it  was  triggered  by  the  first  BSA.  Both  the  pedestal  and  envelope 
of  the  Doppler  signal  were  used  for  burst  detection.  The  FiFo  buffer  mode,  referred  to  “first  in, 
first  out”  mode,  allowed  simultaneous  transmission  and  storage  of  data.  The  FiFo  buffer  was  set 
at  4000  bursts.  Individual  Doppler  bursts  were  counted  because  the  data  collection  mode  was  set 
to  burst.  When  4000  bursts  were  detected,  the  data  was  transferred  to  storage,  and  another  4000 
bursts  were  acquired.  This  process  was  repeated  until  the  maximum  number  of  requested  bursts 
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was  reached,  or  the  timeout  setting  had  expired.  The  quality  factor  was  the  ratio  between  the 
envelope  and  pedestal  levels.  A  Doppler  signal  was  accepted  if  it  had  a  quality  factor  equal  to  or 
greater  than  the  specified  value.  The  oversize  rejection  parameter,  which  rejected  oversized 
particles,  was  disabled  for  this  experiment. 


Table  5-6:  BURSTware  Softkey  Menu  Settings 


BSA#1 

BSA  #2 

Timer  Clock 

Master 

Slave 

Coincidence  Mode 

Master 

Slave 

Arrival  Time  Base 

Internal 

Internal 

Burst  Detection 

Both 

Both 

Oversize  Rejection 

1 

1 

Buffer  Mode 

FiFo 

FiFo 

Max.  Anode  Current 

1.6 

1.6 

Quality  Factor  (%) 

50 

50 

Collection  Mode 

Burst 

Burst 

The  Validation  option  allowed  the  specification  of  what  type  of  data  was  to  be  recorded 
by  each  BSA.  Only  two  options  were  available  for  this  parameter:  All  Data,  or  Valid  Data.  Only 
valid  data  was  recorded  by  each  BSA  for  this  experiment. 

Once  these  parameters  were  entered,  data  collection  could  take  place.  The  high  voltage 
to  the  PM  tubes  was  enabled  at  this  time.  The  Start  option  initiated  data  acquisition  according  to 
the  given  parameter  settings.  The  BURSTware  was  now  ready  to  make  LDV  measurements  with 
the  current  settings. 

5.3.3  BURSTware  Post-processing 

Once  Doppler  burst  data  had  been  collected,  the  Process  option  of  the  software  was  used 
to  initiate  and  monitor  data  processing  and  conversion.  The  type  of  data  to  be  computed  was 
selected  with  the  Setup  option.  For  this  experiment,  velocity  data,  coincidence  filtering,  and 
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moment  data  files  were  processed.  The  coincidence  window  was  set  at  0.000333  msec.  The 
effects  of  the  size  of  the  coincidence  window  were  evaluated  by  varying  it  one  order  of 
magnitude  above  and  below  0.000333  msec,  with  no  apparent  change  in  results.  An  optical 
transform  was  not  required  because  the  optics  were  measuring  the  actual  components  of  velocity. 
After  processing,  both  pre-  and  post-coincidence  moment  files,  which  included  all  mean  and 
turbulent  statistics  from  the  data  collection  run,  were  exported.  The  default  file  extension,  .1st, 
was  used.  These  moment  files  would  be  used  later  by  other  FORTRAN  processing  codes  [23]  to 
determine  turbulent  statistics,  and  flow  properties.  Post-processing  took  about  5  minutes,  which 
provided  enough  time  to  prepare  the  LDV  system  for  data  collection  at  the  next  traverse  station. 
In  other  words,  during  the  time  required  to  empty  the  vacuum  tanks,  the  data  could  be  processed 
and  reviewed  before  the  traverse  was  moved  to  the  next  data  location. 

5.4  Overview  of  a  Typical  Data  Profile 

The  following  section  outlines  the  experimental  procedure  for  data  collection  at  a  single 
traverse  station.  While  waiting  for  the  vacuum  tanks  to  be  evacuated,  the  measurement  volume 
was  moved  to  the  proper  station.  The  wind  tunnel  was  ready  to  run  when  the  vacuum  pressure 
was  below  15  mmHg.  With  the  parameter  file  prepared,  the  computer  was  ready  to  acquire  data. 
The  first  step  in  running  the  tunnel  was  to  open  the  vacuum  gate  valve.  Then,  the  high-pressure 
solenoid  valve  was  opened.  This  order,  first  vacuum,  and  then  high-pressure  was  imperative. 
Once  the  tunnel  was  started,  two  jets  of  the  atomizer  where  opened  introducing  seed  particles 
into  the  flow  field.  A  delay  time  of  about  5  seconds  allowed  the  flow  field  to  stabilize.  The 
computer  data  acquisition  was  begun  after  this  delay  period.  Doppler  bursts  were  collected  for 
12  seconds.  The  minimum  number  of  bursts  collected  was  6000.  However,  in  some  cases,  as 
many  as  135000  bursts  were  collected.  For  most  of  the  data,  approximately  65000  bursts  were 
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collected.  After  the  data  collection  period  expired,  the  atomizer  was  shut  off.  Finally,  the  tunnel 
was  shut  down,  by  closing  the  solenoid  switches  in  reverse  order.  Turning  the  seed  particles  off 
before  shutting  off  the  tunnel  evacuated  all  of  the  olive  oil  from  the  tunnel.  This  was  found  to 
reduce  the  oil  sediments  on  the  windows.  During  the  tunnel  recharge  time,  the  newly  acquired 
data  was  processed. 

This  same  procedure  was  repeated  until  the  boundary  layer  had  been  sufficiently 
mapped.  It  should  be  noted  that  only  two  parameter  files  needed  to  be  created;  one  for  the 
single-component  locations,  and  one  for  the  two-component  stations.  Because  the  first  four  data 
acquisition  locations  of  each  profile  were  in  close  proximity  of  the  wall,  only  single  component 
velocity  measurements  were  collected.  The  green  beams  were  clipped  by  the  wall,  and  an 
accurate  signal  was  not  possible.  Thus,  BSA  #2  was  disabled,  and  the  parameter  files  for  these 
single  component  runs  were  configured  to  only  measure  the  u-component  of  velocity  along  the 
wall.  These  two  files  were  used  as  templates  for  the  remaining  parameter  files.  In  this  manner, 
the  numerous  BSA  settings  were  copied  into  a  new  parameter  file  without  the  need  to  tediously 
re-enter  each  setting. 

The  deposit  of  olive  oil,  and  other  particles,  on  the  glass  windows  interfered  with  the 
laser  beams,  by  causing  the  light  to  be  scattered.  Hence,  the  windows  of  the  tunnel  were 
periodically  removed  and  cleaned  with  distilled  water  and  a  very  soft  cloth,  to  ensure  that  the 
optical  glass  was  not  scratched. 

Vibrations  in  the  floor  due  to  the  vacuum  pumps  caused  slight  changes  in  the  micrometer 
settings;  resulting  in  misalignment  of  the  receiving  and  transmitting  probe  volumes.  Thus,  the 
lasers  were  periodically  realigned  according  to  the  previously  mentioned  procedure  when  this 
occurred. 
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6.  Data  Reduction  Procedures 


6. 1  Calculation  of  Velocity 

The  Doppler  frequency,  as  calculated  by  the  Burst  Spectrum  Analyzer,  is  converted  into 
the  velocity  of  the  particle  by  using  the  frequency- velocity  calibration  factor,  .  The 

calibration  factors,  as  given  in  Chapter  5,  for  the  u-  and  v-component  of  velocity  are 


C.,,.  =8.128-^,  =7.709-"'/^ 


MHz 


MHz 


(6.1) 


Recall,  these  factors  are  based  on  the  optical  setup  and  wavelength  of  the  laser  beams  used. 
Once  the  instantaneous  velocities,  u  and  v ,  are  measured  ,  the  mean,  u  and  v  ,  and  fluctuating, 
u'  and  v' ,  components  of  velocity  may  be  calculated. 


6.2  Calculation  of  Mach  Number,  Temperature,  Density,  and  Pressure 

The  Mach  number  at  each  traverse  station  was  calculated  with  the  local  measured 
velocity  and  total  temperature.  The  flow  was  assumed  adiabatic.  For  turbulent  flow,  mean  and 
fluctuating  terms  are  substituted  into  the  relationship  between  temperature,  stagnation 
temperature,  and  Mach  number.  The  adiabatic  relation 
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becomes 


T„  +  r 
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when  the  expanded  temperatures  and  Mach  number  are  included.  Based  on  the  assumption  of 
small  fluctuating  Mach  numbers  and  the  Strong  Reynolds  Analogy  (SRA),  Morkovin  [28] 
derived  a  first  order  approximation  for  the  ratio  of  fluctuating  and  mean  temperature  as 
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The  Strong  Reynolds  Analogy,  based  upon  the  similarity  between  the  compressible  momentum 
and  energy  equations  when  Pr  =1,  states  that  T'jT^  is  zero  for  adiabatic  flows.  With  these 
assumptions,  the  local  Mach  number  may  be  calculated  by 


M  = 


2Q^ 


1/2 


2YRr,-e“(r-i) 


The  Reynolds  averaged  total  velocity,  Q  ,  is  given  by 


(6.5) 


Q  =  yju^  +y^  (  6.6) 

This  local  Mach  number  is  calculated  assuming  a  constant  total  temperature,  ,  through  the 
boundary  layer.  Kistler  [20]  verified  this  assumption  for  the  flat  plate  case.  He  found  that 
within  the  supersonic  boundary  layer  ( e  1.7,3.56), 


where  the  adiabatic  wall  temperature, 
r  ~  ,  and  the  relationship 
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7^ ,  is  calculated  using  the  adiabatic  recovery  factor. 
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r  = 


T  -T 

aw  e 

T  -T 

^0,~  e 


(6.8) 


with,  ,  being  the  temperature  at  the  edge  of  the  boundary  layer.  Assuming  Pr  =  0.7 1 ,  the 
recovery  factor  is  r  =  0.892  .  Eq.  (6.7)  may  be  rewritten,  in  terms  of  the  Mach  number  at  the 
edge  of  the  boundary  layer,  M  ,  as 


71,  l  +  r[(T-l)/2]M,^ 

1  +  [(Y-1)/2K 


(6.9) 


For  the  present  flat  plate  case,  =  2.9,  and  7^ ^  =  0.932  from  Eq.  (6.9).  Thus, 

~  >  which  is  nearly  constant. 

To  calculate  the  density,  the  assumption  that  the  pressure  gradient  normal  to  the  surface, 
dp/dn ,  equals  zero,  was  used.  While  this  assumption  is  valid  for  the  flat  plate  case,  it  does  not 
necessarily  hold  for  flows  with  adverse  pressure  gradient.  Streamlines  near  a  curved  wall  surface 
will  experience  curvature,  and  a  pressure  gradient  normal  to  the  streamlines  will  be  generated. 
According  to  Spina  [28],  the  magnitude  of  these  gradients  is 


dp  _  —pu^  _  —jM^p 
dn  R  R 


(6.10) 


where  R  is  the  radius  of  curvature.  The  difference  in  pressure  between  the  wall  and  boundary 
layer  edge  is 


(6.11) 


The  change  in  pressure  at  the  x  =  68  cm  APG  station  was  calculated  as  12500  Pa.  At  the 
X  =  71  cm  APG  location,  the  change  in  pressure  across  the  boundary  layer  was  -27500  Pa.  The 
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edge  pressure  was  calculated  by  assuming  that  no  shocks  existed  in  the  freestream.  Thus,  with 
constant  total  pressure  along  the  boundary  layer  edge,  the  edge  pressure  is  given  as 


Pe=Po 


-T 

Y-1 


(6.12) 


where  was  measured  as  described  above. 

By  assuming  constant  total  temperature  and  pressure  across  the  boundary  layer,  the 
Crocco-Busemann  approximation  can  be  used  to  determine  the  mean  temperature  and  density, 
T  and  p  ,  based  on  the  temperature  and  density  at  the  wall,  and  p  ^  [32] : 
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where,  the  coefficients  are  defined  as 
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Using  the  perfect  gas  law  and  the  wall  temperature,  the  density  at  the  wall  is  given  by 


where  R  is  the  gas  constant. 
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(6.16) 
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6.3  Calculation  of  Boundary  Layer,  Displacement,  and  Momentum  Thickness 

Boundary  layer  thicknesses  based  on  both  the  Mach  number  ( 6  ^  ),  and  u-component  of 

velocity  (5„ )  were  calculated.  These  parameters  were  calculated  to  enable  comparison  between 
conventional  mean  flow  probe  (Pitot  and  cone-static  pressure)  data  [9]  which  result  in  direct 
measurements  of  the  Mach  number  and  the  present  LDV  data.  In  this  experiment,  each  boundary 
layer  thickness  was  based  on  99.5%  of  the  edge  value.  Thus,  the  velocity  and  Mach  number 
profiles  were  non-dimensionalized  by  the  edge  values,  and  the  location  above  the  wall  where 
each  quantity  was  99.5%  of  the  edge  value  was  determined.  The  compressible  displacement,  5  * , 
and  momentum,  0  ,  thicknesses  are  defined,  respectively,  as 
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The  incompressible  displacement  thickness,  6*  ,  and  momentum  thickness,  0, ,  were  also 
calculated  by  setting  the  density  ratio,  p/p^ ,  in  Eqs.  (6.17)  and  (6.18)  equal  to  one. 


(6.18) 


6.4  Calculation  of  Turbulent  Statistics 

As  previously  noted  in  Chapter  2,  the  characteristics  of  turbulent  flow  are  analyzed  by 
separating  the  flow  properties  into  a  mean  and  fluctuating  component.  Time-averaged  quantities 
are  used  in  the  turbulence  analysis  because  the  randomness  of  turbulence  makes  exact  calculation 
of  flow  properties  impossible.  By  introducing  the  separated  mean  and  fluctuating  quantities. 
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statistical  correlations  involving  the  fluctuations  and  their  products  are  introduced  into  the 
equations  of  motion  [35],  These  statistical  methods  are  used  to  help  provide  plausible  solutions 
to  the  governing  equations. 

These  turbulent  statistics  also  help  to  characterize  turbulent  flow.  The  fluctuating 
component  may  be  positive  or  negative;  its  time  average  is,  by  definition,  equal  to  zero.  The 
average  size  of  the  fluctuating  term  is  characterized  by  the  root-mean-square  (RMS)  value  of  the 
fluctuation.  The  RMS  value,  or  standard  deviation,  is  defined  as 


(6.19) 


where  there  are  n  samples  in  the  summnation  term.  The  variable,  x ,  may  represent  either  u-  or 
v-component  of  velocity.  The  RMS  value  is  also  used  to  determine  higher-order  moments  such 
as  the  skewness  and  flatness  factors.  The  skewness  ( Sk^ ),  and  flatness  ( Fl^ )  factors  are  given 


by 
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The  skewness  factor  is  a  method  of  quantifying  to  which  side  of  the  mean  the  majority  of  data 


points  are.  Negative  skewness  values  indicate  that  more  values  are  below  the  mean.  These  two 
turbulent  statistics  help  determine  if  there  are  seeding  biasing  problems  in  the  LDV  results  [12]. 


Equations  (6.19)  -  (6.21)  are  not  the  unbiased  statistical  estimates.  The  denominator  should  be 
divided  by  n  —  1  instead  of  n  .  However,  since  the  minimum  number  of  points  exceeded  6000 
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for  all  data,  this  resulted  in  less  than  a  0.01%  bias.  Thus,  it  was  deemed  acceptable  to  use  the 
BURSTware  results. 

The  intermittency  factor,  y  „ ,  is  another  valuable  turbulence  statistic,  that  helps  to 
describe  the  flow.  The  edge  of  a  turbulent  boundary  layer  is  not  always  turbulent.  The 
turbulence  becomes  intermittent,  where  there  are  intervals  of  time  when  the  flow  is  not  turbulent. 
The  intermittency  factor  is  the  ratio  of  the  flattening  factors  in  the  wholly  turbulent  region  to  that 
of  the  intermittent  region.  Near  the  wall,  where  the  flow  is  wholly  turbulent,  the  flattening  factor 
is  near  a  constant  value  of  3.0  [21],  Thus,  the  intermittency  factor  is  given  by 


and 


( 6.24) 
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The  velocity  correlation,  u'v' ,  is  one  of  the  new  terms  that  results  in  the  RANS 
equations.  When  multiplied  by  the  average  density,  this  term  becomes  the  first  term  in  the 
Reynolds  shear  stress,  or  as  discussed  in  Chapter  2  [Eq.  (2.13)],  the  Havre  shear  stress.  The 
velocity  correlation  is  given  by 


u'v'  =  '^{u'v')jn 


The  velocity  correlation  coefficient,  ,  is  defined  as 


( 6.26) 
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This  term  can  be  used  as  a  measure  of  the  anisotropy  of  the  turbulent  flow.  For  isotropic  flow, 
=  1.0 .  The  Cauchy-Schartz  [2]  inequality  dictates  that  |/?„J  ^  1.0.  For  incompressible 


boundary  layers,  this  term  is  nominally  0.5  [3]. 

The  turbulent  kinetic  energy  (TKE)  provides  a  measure  of  the  turbulent  energy  within 
the  boundary  layer.  TKE  is  defined  as 


TKE  = 


2 


( 6.28) 


The  w-component  of  velocity  could  not  be  measured  with  the  configuration  of  the  LDV  for  the 
present  study.  Thus,  the  fluctuating  velocity,  w' ,  could  not  be  measured.  To  estimate  the  value 
of  the  TKE,  it  was  assumed  that  the  fluctuations,  v'  and  w' ,  were  the  same  order  of  magnitude. 
Hence,  the  estimated  TKE  was  given  by 
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Each  of  these  turbulent  quantities  and  statistics  was  calculated  at  each  traverse  station 
across  the  boundary  layer.  The  results,  and  comparisons,  with  the  data  in  the  literature  are  given 
in  Chapter  7. 


6.5  Calculation  of  Wall  Shear  Stress 

The  flat  plate  wall  shear  stress,  T  is  an  important  parameter  that  will  be  used  later  for 
the  inner  region  (law  of  wall)  van  Driest  velocity  profile.  Van  Driest  [3]  gives  the  relation  as 


0.242 


Sin 


+  sin 


-1 


B' 


V4A'^+fi'^ 


( 6.30) 
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where  K  =  0.41 ,  CO  =  0.68 ,  and  A'^  and  B'  are  as  previously  defined.  Eq.  (6.30)  must  be 
solved  iteratively  for  the  skin  friction,  .  Once  the  skin  friction  has  been  calculated,  the  wall 

shear  stress  may  be  calculated  using 
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6.6  Calculation  of  Van  Driest  Velocity  Profile 

Van  Driest  [32]  proposed  a  transformation  to  the  velocity  profile,  so  that  compressible 
flows  could  be  compared.  The  effective  velocity,  ,  is  given  by 
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( 6.32) 


The  effective  velocity  is  non-dimensionalized  by  the  wall-friction  velocity,  u* ,  where 
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u  = 


( 6.33) 


The  non-dimensional  velocity,  Ju* ,  is  plotted  versus  another  scaled  y-coordinate  given  as 


( 6.34) 


where  is  the  fluid  viscosity  calculated  at  the  wall  temperature.  White  [34]  gives  an  empirical 
relationship  between  these  two  quantities  for  the  flat  plate  case: 
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( 6.35) 


with  K  =  0.41  and  C  =  5.0 .  Cole’s  wake  parameter,  IT  ,  is  given  by 


n  =o.8(p+ojf'' 

where  p  is  Clauser’s  equilibrium  parameter. 


( 6.36) 


P«— ^  (6.37: 

T,  dx 

Equation  (6.35)  incorporates  the  overlap  log  law  formulation  with  Coles  law  of  the  wake 
term.  This  equation  was  developed  for  low-speed  data,  but  the  van  Driest  transformation 
collapses  high-speed  data  to  this  form  as  well.  The  region  below  ~  10 ,  where  the  flow 
becomes  laminar,  is  called  the  laminar  sublayer.  In  this  region,  u'^  =  y'^ . 
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6.7  Calculation  of  Strain  Rates 

Central-difference  strain  rates  were  calculated  on  the  APG  model.  The  central  location 
(i,j),  was  located  71  cm  from  the  throat  of  the  nozzle.  Figure  6-1  shows  the  finite-difference  grid. 


y 


Figure  6-1;  Finite-Difference  Grid 
The  central-derivatives  at  the  point  are  estimated  as  [1]: 


du 

dy 


i.i.k 


yi,j+l,k  yi,J-\,k 


( 6.38) 
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( 6.39) 
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( 
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1 

( 6.40) 


where  dx^  =  J  and  dx^  =  —  ^ij.k  ■  The  strain  rates  for  the  v-component  of 

velocity  may  also  be  calculated  by  substituting  v  for  u  in  Eq.  (6.38)-(6.40). 
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7.  Experimental  Results 


The  results  of  this  investigation  are  presented  in  this  chapter.  Laser  Doppler 
Velocimetry  measurements  were  taken  for  both  the  adverse  pressure  gradient,  using  the 
compression  ramp  model,  and  the  zero  pressure  gradient,  or  flat  plate  case.  These  measurements 
were  made  at  three  locations  downstream  of  the  nozzle;  at  x  =  44  and  7 1  cm  for  the  flat  plate, 
and  at  X  =  68  and  71  cm  for  the  compression  ramp  model.  Both  the  flat  plate  mean  flow  and 
turbulence  results  of  the  present  study  are  compared  to  those  results  of  other  researchers.  In  this 
way,  the  LDV  methods  used  in  the  present  study  may  be  validated.  Then,  the  mean  flow  and 
turbulence  results  of  the  compression  ramp  model  are  compared  to  the  flat  plate  results  to 
determine  the  effect  of  the  adverse  pressure  gradient  on  the  turbulent  boundary  layer.  Table  7-1 
summarizes  the  boundary  layer  edge  conditions  for  all  three  measurement  locations  and  test 
models. 


Table  7-1:  Boundary  Layer  Edge  Conditions 


Flat  Plate 

X  =  44  cm 

Flat  Plate 

X  =  71  cm 

Adverse  Pressure 

X  =  68  cm 

Favorable  Pressure 

X  =  71  cm 

2.72 

2.51 

596.0  m/s 

575.0  m/ s 

Te 

lll.OK 

1 16.0  K 

120.0  K 

BLOK 

Pe 

6823.0  Pa 

7866.0  Pa 

8917.0  Pa 

11985.0  Pa 

P. 

0.213  kg  jm^ 

0.237  kglrn^ 

0.258  kg/rn^ 

0.319%/m^ 

dp^ldx 

0 

0 

1250  Pa/ cm 

-1900  Pa/ cm 

P 

0 

0 

1.12 

-0.94 

7. 1  Flat  Plate  Results 

LDV  measurements  were  made  on  the  flat  plate  model  in  the  same  location  as  those 
made  in  the  pressure  gradient  case  for  two  reasons.  The  first  reason  was  to  provide  a  method  of 
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comparing  the  flat  plate  LDV  results  from  the  present  study  to  those  of  other  accepted  data  from 
the  literature.  A  few  of  these  studies  have  incorporated  LDV  as  the  measurement  technique  (Ref. 
[12,  18-19,  25]).  The  second  reason  for  making  LDV  measurements  on  the  flat  plate  at  the  same 
location  was  to  provide  a  means  of  assessing  the  effects  of  the  adverse  pressure  gradient  on  the 
turbulent  boundary  layer  from  those  of  the  flat  plate  at  the  same  location.  In  this  manner,  the 
effects  of  the  pressure  gradient  could  be  distinguished  from  the  flat  plate.  Using  the  same  AFIT 
facilities  and  compression  ramp  model,  Dotter  [9]  used  Pitot  and  cone  static  methods,  and  hot¬ 
wire  anemometry  to  perform  a  similar  investigation.  Hence,  the  present  data  will  be  compared  to 
those  of  Dotter  as  well. 

7.1.1  Flat  Plate  Mean  Flow  Results 

The  flat  plate  LDV  profile  measurements  were  made  at  44  and  7 1  cm  downstream  of  the 
nozzle.  Both  of  these  locations  were  chosen  to  compare  with  the  flat  plate  and  APG  results  of 
Dotter.  The  edge  velocities  at  the  44  and  71  cm  flat  plate  locations  were  605.0m/5  and 
602.0 m/5,  respectively.  The  velocity  boundary  layer  thicknesses  based  on  a  99.5%  criteria  were 
5.73  mm  and  9.86  mm. 
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Figure  7-1;  Flat  Plate  Velocity  Profile,  x  =  44  and  71  cm 

The  velocity  profiles  were  non-dimensionalized  by  the  freestream  velocity.  Figure  7-2 
presents  the  non-dimensionalized  velocity  profiles  versus  y/5  .  Also  shown  are  the  Pitot-cone 
static  results  from  Dotter  [9].  As  can  be  seen,  the  agreement  is  excellent.  Note  that  Dotter’s 
results  were  scaled  by  6  ^  ,  the  Mach  number  boundary  layer  thickness.  In  order  to  correctly 
compare,  those  results  were  rescaled  to  the  present  boundary  layer  thickness. 

Unfortunately,  LDV  measurements  with  the  optical  quality  windows  at  44  cm  were  not 
possible  because  the  windows  could  not  be  placed  at  this  location.  Hence,  the  measurements 
were  taken  through  1.91  cm  thick  Plexiglas.  It  was  anticipated  that  the  LDV  results  would  be 
distorted.  To  the  contrary,  the  LDV  results  matched  Dotter’s  conventional  pressure  probe 
results,  as  is  evidenced  by  the  profiles  in  Figure  7-2. 


7-3 


Figure  7-2:  Flat  Plate  Non-dimensional  Velocity  Profiles,  x  =  44  and  71  cm 

The  flat  plate  Mach  number  profiles  at  both  locations  are  presented  in  Figure  7-3.  At  44 
cm,  the  edge  Mach  number  measured  2.9  and  the  Mach  based  boundary  layer  thickness,  5  ^  ,  was 
8.20  mm.  At  71  cm,  the  edge  Mach  number  was  2.8,  and  5  ^  measured  10.82  mm.  As 
expected,  the  boundary  layer  thickness  over  the  flat  plate  increased  between  the  locations.  The 
difference  in  Mach  number  at  the  two  stations  was  the  result  of  increased  boundary  layer 
thickness  in  the  constant  area  test  section. 
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Figure  7-3:  Flat  Plate  Mach  Number  Profiles,  x  =  44  and  71  cm 

7.1.2  Flat  Plate  Turbulence  Results 

7. 1.2.1  Flat  Plate  Turbulence  Intensities 

The  u-  and  v-velocity  turbulence  intensities  are  presented  in  Figures  7-4  and  7-5.  Elena 
and  LaCharme  [12],  and  Robinson  et  al.  [25]  made  LDV  measurements  in  Mach  2.9  flat  plate 
flow.  Kistler  [20]  made  hot-wire  measurements  in  the  listed  range  of  Mach  numbers.  The 
Klebanoff  curve  [21]  shows  the  incompressible  turbulence  intensities  measured  with  hot-wire 
anemometry  for  flat  plate  flow.  Finally,  Dotter’s  [9]  x  =  44  cm  cross-wire  results  are  included. 

In  particular,  comparisons  with  the  results  of  Elena,  Robinson,  and  the  present  study  are 
important  because  the  LDV  technique  was  used  in  each  case.  For  both  turbulence  intensities,  the 
agreement  with  the  present  results  is  remarkable.  The  present  LDV  results  very  closely  resemble 
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the  results  of  Elena  and  LaCharme.  Johnson’s  v-component  turbulence  intensity  results  match 
slightly  better  with  the  present  data  than  the  u-component  results.  The  Klebanoff  hot-wire  data 
also  agrees  very  well  with  the  present  data.  However,  the  Kistler  and  Dotter  turbulence  intensity 
levels  were  consistently  off.  This  could  be  a  result  of  the  convenient  assumption  used  in  both 
studies  that  the  static  pressure  fluctuations  are  zero.  These  figures  also  indicate  that  the 
freestream  turbulence  of  the  AFIT  Mach  3.0  wind  tunnel  was  approximately  1.0%. 


Figure  7-4:  Flat  Plate  u-component  Turbulence  Intensities 
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Figure  7-5;  Flat  Plate  v-component  Turbulence  Intensities 


7. 1.2.2  Flat  Plate  Reynolds  Shear  Stress 

The  turbulent  shear  stresses  for  the  flat  plate  at  71  cm  are  presented  in  Figure  7-6.  Note, 
to  scale  the  present  results  as  shown  in  Figure  7-6,  the  van  Driest  skin  friction  law  and  Crocco 
integrals  (see  Chapter  6)  were  used.  The  wall  shear  stress  was  estimated  as  approximately  67  Pa 
for  both  flat  plate  measurement  locations.  As  can  be  seen,  the  turbulent  shear  stresses  in  the  flat 
plate  tend  towards  the  wall  shear  stress  value,  t  ,  as  expected.  The  present  results  agree  very 
well  with  the  cross-wire  results  of  Dotter.  The  difference  between  the  two  results  are  well  within 
the  expected  experimental  scatter. 

In  addition,  the  present  data  also  agrees,  at  least  qualitatively,  with  the  results  of  other 
researchers.  The  flat  plate  LDV  results  of  the  other  researchers  all  exhibit  a  “roll  off’  in  the 
turbulent  shear  stress  measurements  near  the  wall.  This  result  is  attributed  to  angular  bias  errors 
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in  the  LDV  measurements  [12],  and  stationary  fringes  in  the  probe  volume  [17],  As  mentioned 

in  Chapter  3,  angular  bias  can  lead  to  low  estimates  of  u'v' .  Thus,  the  low  fluctuating  velocity 
cross-correlation  causes  low  estimates  in  the  turbulent  shear  stress.  Because  the  optics  were 
aligned  to  avoid  the  problems  of  angular  bias,  and  a  Bragg  cell  was  used  to  generate  moving 
fringes,  the  turbulent  shear  stress  values  for  the  present  flat  plate  LDV  measurements  approach 
the  predicted  wall  shear  stress  values. 


-pu'v' 

Figure  7-6:  Flat  Plate  Turbulent  Shear  Stress 


The  intermittency  function,  which  is  a  measure  of  the  percentage  of  time  that  the 
boundary  layer  is  turbulent  (see  Chapter  6),  is  presented  in  Figure  7-7.  By  generating  a  curve  fit, 
the  constants  used  in  the  intermittency  equation  (Eq.  6.24),  and  thus  the  mean  boundary  layer 
position,  about  which  the  turbulence  fluctuates,  may  be  determined.  For  the  flat  plate  case  at  x  = 
71  cm,  the  standard  deviation  was  0.25  ,  or  1.94  mm  about  the  mean  location  of  y/6  =  0.62  . 
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Thus,  the  flow  is  intermittently  turbulent  at  62%  of  the  boundary  layer  thickness.  The  results 
from  the  present  study  resemble  those  of  Elena  and  LaCharme  [12]  quite  well.  The  scatter  in  the 
present  data  above  y/5  =  1.0  may  be  attributed  to  the  freestream  tunnel  turbulence.  Near  the 
wall,  the  flatness,  Fl^  ,  increased  due  to  the  wide  range  of  measured  velocities. 


Figure  7-7:  Flat  Plate  Intermittency  Function 

In  summary,  this  section  demonstrates  that  the  LDV  technique  used  in  the  present  study 
was  providing  consistent  and  accurate  results.  The  Mach  number  profile,  turbulence  intensities, 
turbulent  shear  stress,  and  intermittency  function  measured  by  the  present  LDV  study  match 
other  accepted  data  quite  well.  The  problems  associated  with  angular  bias  were  avoided  by  the 
optical  arrangement  used  in  the  present  study.  Table  7-2  summarizes  the  computed  flat  plate 
boundary  layer,  displacement,  and  momentum  thicknesses.  With  the  flat  plate  comparisons 
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shown,  the  next  section  discusses  the  effects  of  the  pressure  gradient  on  the  turbulent  boundary 
layer.  These  pressure  gradient  results  were  compared  to  the  flat  plate  results. 


Table  7-2:  Flat  Plate  Boundary  Layer  Thickness 


Location 

8„  (mm) 

5’  (mm) 

0  (mm) 

bbbi 

X  =  44  cm 

5.74 

8.20 

0.80 

0.920 

0.352 

0.369 

X  =  71  cm 

9.86 

10.82 

1.31 

1.27 

0.500 

0.615 

7.2  Pressure  Gradient  Results 

Laser  Doppler  Velocimetry  measurements  were  made  at  two  locations  on  the 
compression  ramp  model,  7 1  cm  and  68  cm.  The  first  location  was  chosen  so  as  to  provide 
comparison  to  Dotter’s  results  at  the  same  location.  By  observing  the  surface  pressure 
distribution  along  the  compression  ramp  model,  it  was  determined  that  even  though  the  x  =  7 1 
cm  location  was  on  the  upstream  side  of  a  compression  ramp,  it  was  actually  in  a  region  of 
favorable  pressure  gradient  (FPG).  Thus,  to  enable  evaluation  of  adverse  pressure  gradient 
effects,  that  were  not  combined  with  the  favorable  pressure  gradient  region,  LDV  measurements 
were  made  3  cm  upstream  at  the  68  cm  location.  Both  sets  of  results  were  also  compared  to  the 
flat  plate  LDV  results  obtained  at  7 1  cm. 


7.2.1  Flow  Visualization 

The  shadowgraph  and  schlieren  images  of  the  compression  ramp  region  are  given  in 
Figures  (7-8)  through  (7-10).  Due  to  the  short  duration  of  the  flash  (10  nanoseconds),  the  flow 
appears  “frozen”,  and  important  qualitative  features  of  the  flow  field  can  be  observed  from  these 
pictures.  The  flow  is  from  left  to  right  in  each  photograph.  Coalescing  compression  waves  can 
be  seen  in  the  flow  over  the  compression  ramp.  These  waves  have  coalesced  into  a  shock  wave 
outside  of  the  boundary  layer.  The  weak  shock  waves  at  the  front  of  the  test  section  were 
produced  by  slight  misalignments  in  the  floors  of  the  test  section.  In  the  schlieren  photographs. 
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the  turbulent  eddies  of  the  attached  boundary  layer  can  clearly  be  seen,  as  well  as  the  non- 
uniform  edge  of  the  turbulent  boundary  layer.  The  thickness  of  the  boundary  layer  decreases  in 
the  streamwise  direction.  From  the  shadowgraph,  the  boundary  layer  thickness  at  x  =  68  and  71 
cm  were  estimated  as  7.52  mm  and  8.67  mm,  respectively.  In  the  schlieren  photographs,  these 
same  thicknesses  were  estimated  at  7.37  mm  and  9.22  mm,  respectively.  Referring  to  Figure  7- 
10,  there  appear  to  have  been  present  very  organized  turbulent  structures,  which  were  aligned  at 
an  angle  of  roughly  40-60°  relative  to  the  wall. 

7.2.2  Wall  Pressure  Distribution  and  Temperature 

Pressure  measurements  were  made  along  the  compression  ramp  model  centerline.  These 
results,  presented  in  Figure  7-11,  show  the  increase  in  surface  pressure  from  the  constant  flat 
plate  value.  The  locations  where  LDV  measurements  were  made  are  indicated  by  the  dashed 
lines.  The  surface  pressure  begins  to  increase  approximately  3  cm  from  the  front  of  the  model. 
Figure  7-11  depicts  the  favorable  pressure  gradient  (FPG)  in  the  vicinity  of  x  =  71  cm.  At  this 
location,  dpidx  <  0 ,  and  the  flow  is  expanding  around  the  top  of  the  model.  The  wall  pressure 
gradient,  dp^jdx ,  at  x  =  68  cm  was  estimated  to  be  1250  Pa/cm.  A  numerical  analysis  of  the 
present  study  [14]  indicated  that  dp^jdx  ~  dp^jdx  for  the  present  flow  field.  These  values 
were  used  in  the  estimation  of  Clauser’s  equilibrium  parameter,  P  ,  listed  in  Table  7-1. 


7-11 


Figure  7-10:  Vertical  Knife  Edge  Schlieren  Photograph 


7-12 


55.0000  60.0000  65.0000  70.0000  75.0000  80.0000 

Distance  from  nozzle  throat,  cm 

Figure  7-11;  Compression  Ramp  Surface  Pressure  Distribution 


The  wall  temperature  measured  295  K.  By  using  Eq.  (6.7)  and  the  freestream  reference 
temperature,  ^  ,  of  298  K,  the  adiabatic  wall  temperature,  ,  was  estimated  as  278  K.  Even 

though  the  calculated  adiabatic  wall  temperature  was  approximately  6%  below  the  measured  wall 
temperature,  the  assumption  of  an  adiabatic  wall  was  used  in  this  analysis. 

7.2.3  Pressure  Gradient  vs.  Flat  Plate  Mean  Flow  Results 

The  mean  flow  results  of  the  LDV  measurements  at  both  x  =  44  and  x  =  71  cm  are 
presented  in  Figures  7-12  through  7-16.  Recall,  the  main  objective  is  to  determine  the  effects  of 
an  adverse  pressure  gradient  on  supersonic  turbulent  boundary  layers.  To  enable  these  pressure 
gradient  effects  to  be  distinguished  from  the  effects  due  to  the  natural  evolution  of  the  flat  plate 
boundary  layer,  both  flat  plate  and  pressure  gradient  results  are  presented  together. 
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The  LDV  measurements  are  compared  with  the  conventional  and  hot-wire  anemometry 
results  of  Dotter  [9],  for  both  the  adverse  and  favorable  pressure  gradient  regions.  Even  though 
Dotter’s  measurements  were  made  in  the  FPG  region,  they  are  still  annotated  as  APG.  Since 
Dotter  used  the  Mach  number  boundary  layer  thickness  to  scale  his  results,  a  scale  factor,  equal 
to  the  ratio  of  the  Mach-  and  velocity-based  boundary  layer  thicknesses,  was  used  to  correlate  the 
data  by  multiplying  it  by  Dotter’s  tabulated  y/6  values.  At  the  71  cm  station,  this  ratio  was 
1.30.  In  all  cases,  the  LDV  mean  flow  results  closely  resemble  those  obtained  in  the  previous 
experiment. 

The  velocity  profiles  are  presented  in  Figures  7-12  and  7-13.  These  figures  show  the 
progression  of  the  velocity  profiles  as  moved  further  down  the  tunnel.  Over  the  compression 
ramp,  the  edge  velocity  decreases  from  600  m/5  to  about  575  m/5  in  a  distance  of  3  cm. 
Meanwhile,  the  velocity  profiles  at  68  cm  and  71  cm,  for  the  APG  and  flat  plate  cases, 
respectively,  are  almost  identical.  In  the  latter  figure,  the  velocity  has  been  non-dimensionalized 
by  the  boundary  layer  edge  velocity.  As  expected,  all  of  the  profiles  collapsed  rather  well. 
Dotter’s  data  at  both  measurement  locations  seems  to  bound  the  LDV  results. 
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The  flat  plate  and  pressure  gradient  density  profiles  are  presented  in  Figure  7-14.  This 
figure  demonstrates  the  excellent  agreement  between  the  LDV  and  conventional  pressure-probe 
results  of  Dotter.  The  compression  ramp  increased  the  density  at  the  boundary  layer  edge.  This 
trend  is  evident  as  the  profiles  decrease  in  magnitude  as  the  distance  down  the  tunnel  is 
increased.  The  greatest  variation  in  profile  shape  occurred  at  x  =  68  cm,  where  the  maximum 
value  is  at  y/b  =  0.6  ,  instead  of  at  the  boundary  layer  edge. 
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Figure  7-14:  Pressure  Gradient  vs.  Flat  Plate  Density 

As  with  the  streamwise  velocity  comparison  above,  which  displays  the  Mach  number 
profiles  for  both  the  flat  plate  and  APG  test  sections  at  each  measurement  location.  Figure  7-15 
shows  the  progression  of  the  boundary  layer.  The  present  results  and  Dotter’ s  results  compare 
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quite  well;  with  the  Mach  number  profiles  at  x  =  44  cm  exhibiting  remarkable  agreement.  The 
flow  slows  from  Mach  2.9  upstream  to  Mach  2.5  at  71  cm  in  the  compression  ramp  region. 
Between  the  measurement  stations  at  x  =  68  and  71  cm,  the  flow  slows  from  Mach  2.72  to  Mach 


1.00  1.50  2.00  2.50  3.00 

Mach  Number 


Figure  7-15:  Pressure  Gradient  vs.  Flat  Plate  Mach  Number 

The  van  Driest  velocity  transformation  helps  to  correlate  compressible  flow  data.  These 
results  were  calculated  and  are  presented  in  Figure  7-16.  The  dashed  line  represents  the  laminar 
sublayer,  where  j u  =  .  The  equation  of  the  solid  line  is  given  as 


* 


u 


2.51n(y^)-h4.9-i-5.12sin' 


(7.1) 
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The  single  component  measurements  were  made  closer  to  the  wall  than  the  regular  two- 
component  LDV  measurements.  With  single  component  LDV,  the  closest  station  to  the  wall  was 
=45.  Overall,  the  measured  LDV  velocity  profile  compared  very  well  with  the  theoretical 
values  predicted  by  van  Driest.  Note,  the  large  uncertainty  is  due  to  the  estimate  of  skin  friction 
(±  10  %,  Ref.  [3]),  and  possible  error  in  the  wall  shear  stress  estimation. 


Figure  7-16:  Van  Driest  Velocity  Profile 

As  with  the  flat  plate,  the  boundary  layer,  displacement,  and  momentum  thicknesses 
were  calculated  in  the  adverse  pressure  gradient  region.  These  results  are  presented  in  Table  7-3. 
When  compared  to  the  flat  plate  values,  it  is  evident  that  the  velocity  and  Mach  boundary  layer 
thicknesses  decrease  as  a  result  of  the  adverse  pressure  gradient.  Both  boundary  layer 
thicknesses  at  the  x  =  71  cm  station  are  smaller  than  those  at  the  68  cm  location.  This  indicates 
that  the  streamlines  are  being  compressed  towards  the  wall,  and  the  APG  effects  are  still  present 
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even  in  the  FPG  region  at  x  =  71  cm.  These  results  agree  well  with  the  analysis  of  Sturek  [30], 
who  showed  that  even  though  the  boundary  layer  thickness  is  decreasing,  mass  is  being  entrained 
in  the  boundary  layer. 


Table  7-3:  Pressure  Gradient  Boundary  Layer  Thickness 


Location 

6„  (mm) 

5‘  (mm) 

0  (mm) 

besbh 

X  =  68  cm 

9.09 

10.04 

1.97 

1.26 

0.672 

0.753 

X  =  71  cm 

7.48 

8.44 

1.95 

0.97 

0.447 

0.565 

7.2.4  Pressure  Gradient  vs.  Flat  Plate  Turbulence  Results 

The  remainder  of  this  chapter  presents  the  comparison  of  the  flat  plate  and  adverse 
pressure  gradient  turbulence  terms.  The  turbulence  intensities,  Reynolds  shear  stresses,  and 
intermittency  function  are  presented,  along  with  the  higher-order  moments,  such  as  flatness  and 
skewness.  Lastly,  the  turbulent  kinetic  energy  (TKE)  and  velocity-velocity  correlation 
coefficient  are  examined.  In  all  cases,  the  flat  plate  results  at  x  =  71  cm,  and  the  results  for  the 
compression  ramp  measurement  locations  are  given  in  the  same  figure  for  ease  of  comparison. 

7. 2.4.1  Turbulence  Intensities 

The  turbulence  intensities,  which  give  a  measure  of  the  magnitude  of  the  fluctuating 
component,  are  presented  in  Figure  7-17  and  Figure  7-18.  The  turbulence  intensities  (T.I.)  for 
both  the  u-  and  v-component  increased  as  a  result  of  the  pressure  gradient,  with  the  streamwise 
component  exhibiting  the  largest  increase.  For  the  u-component,  the  maximum  turbulence 
intensity  at  x  =  68  cm,  located  at  y/6  =  0.25 ,  was  9.0%.  Thus,  the  adverse  pressure  gradient 
caused  the  streamwise  turbulence  intensity  to  increase  by  24%  of  its  flat  plate  value,  at  the  same 
location  in  the  boundary  layer.  For  the  FPG  location,  the  u-component  turbulence  intensity 
showed  a  14%  increase  over  the  flat  plat  value.  For  the  v-component,  the  T.I.  at  both  68  and  71 
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cm  have  a  maximum  value  of  3.7%  at  y/b  =  0.2 ,  which  results  in  an  8.1%  increase  over  the 
corresponding  flat  plate  value.  Thus,  the  FPG  effects  have  reduced  the  magnification  of  the 
turbulence  intensities. 

It  should  be  noted  that  the  results  of  the  hot-wire  analysis  were  very  different.  Above 
y/b  ~  0.5 ,  the  hot-wire  and  LDV  u-component  T.I.  results  match  quite  well.  Below  this 
location,  the  hot-wire  results  are  lower  than  those  of  the  present  study.  From  Figure  7-18,  it 
appears  as  if  the  hot-wire  turbulence  intensity  exhibits  the  same  general  trend  as  the  LDV  results. 
Yet,  the  hot-wire  results  are  magnified  by  as  much  as  163%  of  the  LDV  values.  These  trends 
seem  to  indicate  that  the  turbulence  transformation  [9]  used  to  extract  the  turbulence  intensities 
from  the  hot-wire  analysis  may  not  be  valid  for  flow  fields  inyolving  pressure  gradients. 
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Figure  7-17;  Pressure  Gradient  vs.  Flat  Plate  u-component  Turbulence  Intensities 
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Figure  7-18:  Pressure  Gradient  vs.  Flat  Plate  v-component  Turbulence  Intensities 
7.2.4.2  Reynolds  Shear  Stresses 

The  Reynolds  shear  stresses,  which  were  extremely  sensitive  to  the  pressure  gradients, 
are  presented  in  Figures  7-19  and  7-20.  In  the  first  figure,  the  Reynolds  shear  stress,  calculated 

as  the  product  of  the  mean  density  and  fluctuating  velocity  cross  correlation,  m'v'  ,  at  each 
traverse  station,  was  non-dimensionalized  by  the  flat  plate  wall  shear  stress,  T  .  Therefore,  the 
71  cm  flat  plate  results  approach  unity  at  the  wall.  At  71  cm,  the  Reynolds  shear  stress  in  the 
FPG  region  increases  by  67%  of  the  flat  plate  value.  The  68  cm  APG  Reynolds  shear  stress 
reaches  a  maximum  value  of  2.6  at  y/5  =  0.25 ,  which  is  an  190%  increase  over  the  flat  plate 
value  of  0.9  at  the  same  location  in  the  boundary  layer.  These  results  agree  with  Ref.  [22]  which 
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showed  the  turbulent  shear  stresses  had  maximum  values  over  the  wall  shear  stress  near 
>;/5  =0.3. 

In  Figure  7-20,  the  Reynolds  shear  stresses  are  non-dimensionalized  by  the  mean 
velocity,  u  .  The  71  cm  APG  profile  seems  to  exhibit  an  inflection  point  near  y/5  =  0.7 ,  which 
may  be  attributed  to  experimental  scatter.  The  profiles  for  the  flat  plate  and  68  cm  APG  case 
have  the  same  slope,  with  the  APG  Reynolds  shear  stress  increased  to  133%  of  the  flat  plate 
value  at  y/6  =  0.2  .  As  expected,  the  Reynolds  shear  stresses  approach  zero  for  locations 
y/6  >  1.0. 

As  with  the  hot-wire  turbulence  intensities,  the  turbulent  shear  stresses  differ  greatly 
from  the  LDV  results.  Because  the  LDV  technique  directly  measures  the  mean,  fluctuating,  and 
cross-correlation  velocity  components,  which  are  used  in  the  estimates  of  the  turbulence 
intensities,  and  the  Reynolds  shear  stresses,  these  values  are  assumed  correct.  No  assumptions  or 
transformations  were  required  to  extract  these  terms  as  with  the  hot-wire  analysis.  These 
variations  indicate  that  the  assumptions  used  in  the  hot-wire  turbulence  transformation  may  not 
be  valid  for  pressure  gradient  supersonic  flow. 
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Figure  7-19:  Reynolds  Shear  Stress 
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Figure  7-20:  Reynolds  Shear  Stress 


7.2.4.3  Skin  Friction  Estimates 


In  the  logarithmic  or  overlap  region  of  the  boundary  layer,  ~'^ly+~y  ■  The 

dx 

results  in  Figure  7-16  suggest  that  the  near  wall  measurements  were  in  the  log  region  (i.e. 

<  300).  The  wall  shear  stress,  X  ^ ,  may  be  estimated  by  assuming 

(7.2) 

Hence,  the  skin  friction  coefficient,  ,  may  be  accurately  estimated  as 
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Table  7-4  summarizes  the  results  of  this  skin  friction  analysis.  The  adverse  pressure 
gradient  doubles  the  skin  friction  at  the  same  location.  Interestingly,  the  skin  friction  of  the  flat 
plate  and  the  FPG  region  are  almost  the  same.  The  skin  friction,  then,  is  reduced  by  50% 
between  the  x  =  68  cm  and  71  cm  compression  ramp  locations.  The  van  Driest  estimate 
predicted  skin  friction  values  larger  than  the  current  estimates.  This  may  be  a  result  of 
overpredicting  the  wall  shear  stress  by  the  methods  described  in  Chapter  6.  These  results  agree 
well  with  those  of  Smith  [27]  and  Spina  [28]. 


Table  7-4:  Skin  Friction  Estimates 


Location 

Flat  Plate,  x  =  7 1  cm 

APG,  X  =  68  cm 

FPG,  X  =  7 1  cm 

Cj ,  current  estimate 

0.00 11 

0.0022 

0.0010 

Cf ,  van  Driest 

0.0016 

0.0016 

0.0016 

7. 2. 4. 4  Intermittency  Function 

The  intermittency  function  at  both  68  cm  and  71  cm  is  presented  in  Figure  7-21.  Recall, 
this  is  a  measure  of  the  time  that  the  boundary  layer  edge  is  turbulent.  Curve  fits  have  been 
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generated  to  determine  the  constants  from  Eq.  (6.24),  which  give  the  standard  deviation  about  the 
mean  boundary  layer  position.  For  the  68  cm  case,  the  standard  deviation  was  0.175  ,  or  1.54 
mm,  about  the  mean  location  of  y/6  =  0.55.  For  the  71  cm  case,  the  standard  deviation  was 
0.266  ,  or  1.94  mm,  about  the  mean  location  of  y/6  =  0.62  .  Thus,  the  mean  turbulent  boundary 
layer  thickness  is  smaller  in  a  region  of  greater  adverse  pressure  gradient.  When  the  71cm 
intermittency  function  is  compared  to  the  flat  plate  boundary  layer  at  the  same  location,  it  is 
interesting  to  note  that  the  standard  deviation  and  mean  location  are  identical.  Also,  the  mean 
boundary  layer  location  at  the  x  =  68  cm  station  is  smaller  than  the  flat  plate  mean  location. 


Figure  7-21:  Pressure  Gradient  Intermittency  Function 
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7.2.4. 5  Flatness  and  Skewness 

These  higher-order  turbulent  statistics  are  presented  in  Figures  7-22  through  7-25.  In  all 
cases,  the  present  flat  plate  values  match  the  Elena  and  LaCharme  [12]  flat  plate  data  quite  well. 
The  flatness  for  both  components  over  the  flat  plate  tend  to  center  around  zero  throughout  the 
boundary  layer.  The  flatness  in  both  the  u-  and  v-components  in  nearly  zero  to  a  location  of 
y/6  =  0.8  .  This  indicates  that  the  flat  plate  LDV  flatness  measurements  are  not  susceptible  to 
much  variation.  But,  these  values  are  magnified  by  the  adverse  pressure  gradient  model; 
especially  at  the  68  cm  location,  where  the  flatness  begins  to  diverge  from  zero  near  y/5  =  0.5 
for  both  components.  The  skewness,  which  helps  to  determine  to  which  side  of  the  mean  the 
majority  of  the  data  lies,  is  negative  for  the  u-component,  yet  positive  for  the  v-component.  In 
each  case,  within  the  boundary  layer,  these  moments  exhibit  the  trend  of  the  flat  plate  values  as 
one  bound,  the  68  cm,  or  greater  adverse  pressure  gradient  as  the  other  bound,  with  the  71  cm 
APG  case  lying  somewhere  in  the  middle. 
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Figure  7-22:  Pressure  Gradient  vs.  Flat  Plate  u-component  Flatness 
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Figure  7-23:  Pressure  Gradient  vs.  Flat  Plate  u-component  Skewness 
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Figure  7-24:  Pressure  Gradient  vs.  Flat  Plate  v-component  Flatness 
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Figure  7-25;  Pressure  Gradient  vs.  Flat  Plate  v-component  Skewness 
7.2.4.6  Velocity-Velocity  Correlation  Coefficient  and  TKE 

These  turbulent  values  are  presented  in  Figures  7-26  and  7-27,  respectively.  Notice  that 
the  data  of  Elena  and  LaCharme  scatters  much  more  than  their  previous  turbulent  results. 
The  velocity-velocity  correlation  coefficient  results  from  the  present  study  exhibit  a  more 
uniform  trend.  The  coefficient  is  increased  to  65%  of  the  flat  plate  value  at  the  location  of 
_y/6  =  0.3  for  the  x  =  68  cm  station.  Yet,  for  the  location  at  71  cm,  the  coefficient  decreases  by 
27%  of  the  flat  plate  value  at  y/6  =  0.2 .  These  results  indicate  that  the  present  flow  field 
exhibits  incompressible  behavior  near  the  wall,  as  — =  0.4  -  0.5  [3]. 

The  TKE,  which  is  a  measure  of  the  turbulent  kinetic  energy  of  the  velocity  fluctuations, 
shows  increased  values  as  the  wall  is  approached.  This  is  expected  as  the  boundary  layer  is  more 
turbulent  in  this  region,  than  near  the  boundary  layer  edge.  The  flat  plate  and  APG  results  show 
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the  same  trend  in  the  slope  at  the  x  =  71  cm  location.  But,  the  gradient  near  the  wall  at  the  68  cm 
station  seems  to  be  larger.  In  each  case,  the  TKE  approach  0.015%  as  predicted  by  Eq.  (6.29), 
and  u-  and  v-component  turbulence  intensities  of  1.0%.  Recall,  that  the  TICE  has  been  calculated 
assuming  the  fluctuations  in  the  w-direction,  although  not  measured  with  the  present  LDV  study, 
are  of  the  same  magnitude  as  the  v-component  fluctuations. 
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Figure  7-26;  Velocity-Velocity  Correlation  Coefficient 
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Figure  7-27:  Pressure  Gradient  vs.  Flat  Plate  TKE 

7.2.5  Pressure  Gradient  vs.  Flat  Plate  Strain  Rates 

Strain  rate  measurements  were  made  in  the  vicinity  of  the  7 1  cm  location  for  both  the 
compression  ramp  and  flat  plate  models.  Five  measurement  profiles  were  required  for  the 
pressure  gradient  station  at  7 1  cm.  With  reference  to  the  finite-difference  grid  in  Figure  6- 1 ,  the 
distances  ric,  and  dXj  were  4.8  mm  and  9.8  mm,  respectively.  The  lateral  measurements  were 
made  5  mm  on  either  side  of  the  tunnel  centerline.  Strain  rates  were  non-dimensionalized  by  the 
edge  velocity  and  boundary  layer  thickness. 

To  determine  the  relative  strength  of  these  strain  rates,  the  measurements  were  quantified 
after  the  analysis  of  Ref.  [27].  The  ratio  of  each  extra  strain  rate,  e ,  to  the  main  strain  rate  in  the 
boundary  layer,  duldy  ,  was  calculated.  The  distortion  is  classified  as  weak  if  the  ratio 
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ej {duj dy)  =  0.01  ,  and  strong  if  ej (dujdy)  ~  0.1 .  For  the  present  study,  these  ratios  revealed 
that  the  interactions  were  strong. 

Figure  7-28  presents  the  rate  of  change  of  the  u-component  of  velocity  with  change  in 
distance  down  the  tunnel.  As  expected,  this  flat  plate  strain  rate  was  calculated  as  zero.  At 
y/5  =  0.1 ,  the  strain  rate  for  the  APG  case  has  increased  to  0.2.  The  largest  gradient  occurs 
below  y/5  =  0.5 . 


-0.1000  0.0000  0.1000  0.2000  0.3000 
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Figure  7-28:  dujdx  Strain  Rate  Comparison 

The  change  in  u-component  over  the  boundary  layer  thickness  is  given  in  Figure  7-29. 
This  figure  indicates  that  this  strain  rate  for  the  pressure  gradient  case  is  almost  the  same  as  the 
flat  plate  case.  The  quantity  du/dy  is  used  in  determining  gradient  transport  formulations  such 
as  the  eddy  viscosity  model  or  Prandtl’s  mixing  length  model  [3]. 
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Figure  7-29:  duldy  Strain  Rate  Comparison 


Prandtl  proposed  that  the  turbulent  fluctuations  could  be  related  to  a  length  scale  and  velocity 
gradient.  The  Prandtl  mixing  length  model  is  given  by 


where  is  the  mixing  length  and  y  is  the  intermittency  function.  The  mixing  length  is  constant 
for  y/5  >  0.22  and  is  given  by  /„  =  0.095  .  For  locations  below  y/5  =  0.22 ,  =  0.4y . 

Using  the  strain  rates  at  71  cm  and  the  mean  flow  results,  the  turbulent  shear  stress  t was 

calculated.  The  results  are  presented  in  Figure  7-30.  The  Prandtl  mixing  model,  which  is 
primarily  an  incompressible  turbulence  model,  grossly  underpredicts  the  turbulent  shear  stress 


du 

dy 


du 

dy 


(7.4) 
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calculated  in  the  present  experiment.  This  demonstrates  the  need  for  valid  data  so  that  accurate 
compressible  turbulence  models  may  be  developed. 
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Figure  7-30:  Simple  Turbulence  Model 


The  strain  rates  for  the  v-component  of  velocity  are  presented  in  Figures  7-31  and  7-32.  The  flat 
plate  streamwise  strain  fate  for  both  the  u-  and  v-components  was  zero,  as  expected.  The  APG 
strain  rate  in  Figure  7-30  exhibits  an  inflection  point  at  y/5  =  0.6 ,  and  increasing  magnitude 
above  the  boundary  layer  edge.  The  flat  plate  strain  rate  in  Figure  7-32  also  has  an  inflection 
point,  but  closer  to  the  boundary  layer  edge.  The  pressure  gradient  increased  the  strain  rates  for 
both  v-component  measurements. 
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Figure  7-31:  dv/dx  Strain  Rate  Comparison 
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Figure  7-32:  dv/dy  Strain  Rate  Comparison 


8.  Conclusions  and  Recommendations 


8. 1  Conclusions 

An  experimental  investigation  on  the  effects  of  adverse  pressure  gradient  on  the  mean  and 
turbulent  flow  structure  of  a  Mach  2.9  boundary  layer  was  performed.  Laser  Doppler  Velocimetry 
was  used  to  measure  the  mean  and  fluctuating  velocity  components,  turbulence  intensities, 
skewness,  flatness,  and  intermittency  in  the  axial  and  transverse  directions. 

The  present  study  was  conducted  in  the  Air  Force  Institute  of  Technology  Mach  3.0  wind 
tunnel.  LDV  profile  measurements  were  made  at  the  same  location  for  a  flat  plate  model,  and 
compression  ramp  model  that  generated  an  adverse  pressure  gradient.  In  this  manner,  the  effects  of 
the  pressure  gradient  could  be  distinguished  from  the  effects  of  the  flat  plate  boundary  layer.  LDV 
profile  measurements  were  also  made  a  few  centimeters  upstream  of  the  location  studied  by  Dotter 
[9].  It  was  determined  that  the  71  cm  measurement  station  was  in  a  region  of  favorable  pressure 
gradient. 

The  flat  plate  LDV  measurements  were  made  for  two  reasons.  The  first  reason  was  to 
provide  enable  comparison  of  the  data  from  the  present  study  with  other  accepted  data  from  the 
literature.  This  comptirison  enabled  validation  of  the  technique  used  in  the  study.  The  second 
reason  was  to  provide  a  means  of  distinguishing  the  effects  of  the  pressure  gradient  on  the 
turbulent  boundary  layer  from  those  of  the  natural  progression  of  the  flat  plate  boundary  layer. 

Rat  plate  measurements  were  taken  44  and  7 1  cm  downstream  of  the  throat  of  the  nozzle. 
The  Mach  number  decreased  from  2.87  at  the  44  cm  location  to  2.79  downstream.  The  boundary 
layer  thickness  increased  from  5.74  mm  to  9.86  mm  at  these  two  stations.  The  flat  plate 
turbulence  results  matched  those  of  previous  researchers  quite  well,  as  evidenced  by  Figures  7-4 
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through  7-7.  The  flat  plate  LDV  results  of  the  present  study  agreed  very  well  with  LDV  results  of 
other  researchers  [12,  19,  25].  The  Reynolds  shear  stress  approached  the  wall  value,  and  did  not 
exhibit  the  “roll  off’  in  the  turbulent  shear  stress,  that  was  experienced  in  the  previous  LDV 
studies.  This  problems  were  avoided  by  aligning  the  optics  with  the  mean  flow,  and  using  a  Bragg 
cell  to  create  moving  fringes  in  the  probe  volume.  At  the  x  =  71  cm  location,  the  boundary  layer 
was  intermittently  turbulent  to  within  about  y/5  =  0.62  . 

Two  Laser  Doppler  Velocimetry  measurements  profiles,  located  68  and  71  cm  downstream 
of  the  throat  of  the  nozzle,  were  taken  over  the  compression  ramp  model.  At  x  =  68  cm,  the  flow 
experienced  an  adverse  pressure  gradient.  At  the  second  location,  the  flow  was  expanding  around 
the  top  of  the  model  to  generate  a  region  of  favorable  pressure  gradient.  The  mean  flow 
compression  ramp  results  indicate  that  the  Mach  number  decreased  from  2.72  at  x  =  68  cm  to  2.51 
at  the  71  cm  station.  The  boundary  layer  thickness  decreased  from  9.09  mm  to  7.48  mm  between 
the  two  locations.  When  compared  to  the  x  =  71  cm  flat  plate  measurements,  the  Mach  number 
and  boundary  layer  thickness  in  the  APG  region  decreased  by  2.5%  and  7.8%,  respectively.  In  the 
FPG  region,  these  quantities  decreased  by  10%  and  24%,  respectively,  from  the  flat  plate  values. 
The  change  in  boundary  layer  thickness  and  flow  structure  was  observed  with  nanosecond 
shadowgraph  and  schlieren  photography. 

The  pressure  gradient  turbulence  results  indicate  that  the  turbulence  intensities  and  shear 
stresses  increased  over  the  flat  plate  values.  The  u-component  turbulence  intensity  in  the  APG 
region  increased  by  24%  of  the  flat  plate  value,  while  in  the  FPG  region  increased  it  increased  by 
14%.  The  v-component  turbulence  intensities  for  both  locations  increased  by  8. 1%  over  the  flat 
plate  value.  The  Reynolds  shear  stresses  increased  by  190%  and  67%  over  the  flat  plate  values  for 
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the  APG  and  FPG  regions,  respectively.  Thus,  the  FPG  reduced  the  amplification  of  the 
turbulence  intensities  and  turbulent  shear  stresses. 

The  skin  friction  estimate  indicated  that  the  skin  friction  was  twice  the  flat  plate  value  at 
the  68  cm  location.  Further  downstream  in  the  FPG  region,  the  skin  friction  was  almost  identical 
to  the  flat  plate  value.  In  the  APG  region,  the  boundary  layer  was  intermittently  turbulent  to  within 
about  y/6  =  055 .  Thus,  the  adverse  pressure  gradient  decreased  the  mean  turbulent  location  as 
compared  to  the  flat  plate  value.  In  the  FPG  region,  the  mean  locations  were  surprisingly  identical. 

Based  on  the  comparison  of  the  hot-wire  anemometry  and  LDV  results,  it  appears  that  the 
turbulence  transformation  [9]  used  to  extract  the  hot-wire  turbulence  data  may  not  be  valid  for 
supersonic  pressure  gradient  flow.  "Oie  hot-wire  and  LDV  results  compare  quite  well  for  flat  plate 
flow.  But,  large  differences  are  evident  for  the  pressure  gradient  region.  The  LDV  technique 
directly  measures  the  mean  and  fluctuating  velocity  components.  Hence,  the  turbulence  intensities 
and  Reynolds  shear  stresses  are  directly  acquired  without  the  need  for  a  transformation  on  the 
measured  data. 

As  mentioned  in  Chapter  1,  this  investigation  was  designed  to  meet  the  criteria  of  Settles 
and  Dodson  [26].  Table  8-1  shows  how  each  of  the  criteria  were  met. 


Table  8.1:  Settles  and  Dodson  Criteria 


Criteria 

Compliance 

Baseline  Applicability 

Mach  2.9  Flow 

Simplicity 

Simple  Geometry  Modeled  in  Ref  [14] 

Specific  Applicability 

Compared  to  CFD  Code  [14]  for  Validation 

Well  Defined  Experimental  Boundary  Conditions 

Tables  7-1  through  7-3,  Appendix  C 

Well  Defined  Error  Bounds 

Appendix  A 

Adequate  Documentation  of  Data 

Appendix  B 

Adequate  Spatial  Resolution 

Figures  in  Chapter  7 
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8.2  Recommendations 


It  is  recommended  that  more  LDV  profile  measurements  are  made  at  locations  between  the 
68  cm  and  71  cm  stations.  In  this  manner,  the  progression  of  the  supersonic  turbulent  boundary 
layer  may  be  more  fully  understood.  Also,  these  results  may  be  compared  with  CFD  results  [14]  at 
the  same  location.  Other  compression  ramp  models,  which  generate  higher  magnitude  adverse 
pressure  gradients,  may  be  studied  and  compared  to  this  model.  Lastly,  an  imposed  adverse 
pressure  gradient  like  the  study  of  Smith  and  Smits  [27]  may  be  evaluated  to  distinguish  the 
effects  of  streamline  curvature  and  the  pressure  gradient. 
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Appendix  A:  Error  Analysis 


This  section  gives  estimates  of  the  uncertainty  of  the  measured  and  calculated  quantities 
of  the  present  study.  This  analysis  follows  the  well-detailed  error  analysis  of  Luker  [23], 

With  any  experiment  it  is  necessary  to  perform  an  error  analysis  to  determine  the  validity 
of  the  results.  The  error  of  a  computed  quantity  depends  on  the  error  of  the  measured  quantities 
that  were  used  in  the  computation.  The  absolute  error.  A/  ,  in  the  quantity  /  is  determined  by 
differentiating  f  hy  n  dependent  variables  x- : 


dx- 


AX; 


(A.l) 


where  /„  is  the  mean  value  of  /  around  which  the  function  is  linearized.  The  error  in  /  may 
be  expressed  as  a  percentage  by  using  =  Af  x  100%.  Thus,  Eqn.  (A.l)  becomes 


1=1 


fo 


fo 


-e. 


(A.2) 


Thus,  e  ^  denotes  the  percent  error  in  the  function  /  .  These  equations  give  the  maximum 

possible  error  assuming  the  errors  are  simultaneously  at  their  maximum  values.  It  is  reasonable 
to  expect  that  the  errors  are  not  all  maximized  at  the  same  time.  Thus,  the  Euclidean  norm  is 
used  to  determine  the  percent  error: 


r 

2" 

_ 1 

a/ 

[ax,. 

f 

J 

(A.3) 


Error  calculations  throughout  this  section  will  be  estimated  by  Eqn.  (A.3)  unless  otherwise 
indicated. 
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The  error  associated  with  the  LDV  velocity,  temperature,  and  pressure  measurements  is 
given  in  the  next  section.  Once  these  errors  are  defined,  their  influence  on  the  calculated  values 
is  estimated. 

A.l  LDV  Velocity  Error 

The  error  in  the  LDV  velocity  measurements  was  attributed  to  four  key  sources.  These 
sources  included  traverse  position  error,  BSA  error,  angular  alignment  error,  and  repeatability 
error.  The  traverse  error  occurred  because  the  BURSTware  didn’t  accurately  record  the  traverse 
position.  For  example,  LDV  measurements  were  made  at  the  location  y  =  0.09,  z  =  -0.76.  When 
moving  the  traverse  to  the  next  profile  station,  the  software  would  read  the  traverse  position  as 
taking  measurements  at  the  location  y  =  0.09,  z  =  -0.75.  Thus,  the  error  in  the  traverse  position 
would  increase  by  0.01  mm  each  time  the  traverse  was  moved.  The  location  error  associated 
with  the  i"'  data  location  in  the  profile,  in  units  of  mm,  was 

Ay,.  =  Ay„  -I-  0.0 l(i  - 1)  (  A.4) 

where  Ay„  is  the  origin  placement  error.  The  origin  placement  error  was  estimated  by 
evaluating  the  possible  uncertainty  in  the  alignment  of  the  optics  with  the  wall.  When  the  lasers 
were  aligned  with  the  wall,  it  was  assumed  that  the  “spot  pattern”  was  placed  within  25%  of  the 
radius  of  the  spot’s  center.  Thus,  the  angular  error.  At]) ,  was  approximately  4.3°  and  the  origin 

placement  error,  Ay„  ,  was  0.074  mm.  This  is  a  rather  large  error  considering  the  rotation  angles 
used  for  LDV  measurements  with  the  compression  ramp  test  section. 

The  BSA  error  in  the  velocity  measurements  resulted  from  limitations  in  calculating  the 
Doppler  frequency.  The  frequency  error  of  the  BSA  is  [7] 


A/o  = 


n{BW) 

2*iV 


(A.5) 


where  N  is  the  record  length,  and  n  ,  based  on  the  bandwidth,  were  1.2  and  1.5  for  BSA  #1  and 
#2,  respectively.  Thus  the  BSA  errors  for  the  u-  and  v-components  of  velocity  were  0.095  m/s 
and  0.090  m/s,  respectively. 

The  repeatability  errors  were  determined  by  making  twenty  test  runs  in  the  wind  tunnel 
(see  Appendix  C).  The  repeatability  errors  are  given  in  Table  A-l._  The  errors  in  the  mean  and 
fluctuating  v-components  of  velocity  were  normalized  by  u  . 


Table  A-1:  Repeatability  Velocity  Errors 


e 

u 

£ 

V 

Ev' 

£  — 

U  V 

0.19% 

0.04  % 

0.07  % 

0.03  % 

0.02  % 

Another  possible  source  of  error  was  the  off-axis  angular  alignment  of  3.5°.  Thus,  the 
lasers  were  measuring  a  lower  velocity  component  than  the  actual  value.  Based  on  the  cosine  of 
the  off-axis  angle,  the  actual  velocity  was  1.002  times  the  measured  value.  The  resulting  error, 
0.2%,  was  ignored  when  compared  to  the  magnitudes  of  the  other  errors. 

The  total  velocity  error  at  each  measurement  location  in  the  profile  is  given  by  the 
Euclidean  norm  of  each  of  those  errors  mentioned  above: 


The  velocity  error  was  calculated  at  each  position  in  the  profile.  To  simplify  further  error 
calculations,  the  velocity  error  was  assumed  constant  and  equal  to  the  maximum  error  in  the 
profile.  Table  A-2  presents  the  LDV  velocity  errors,  where  again  the  v-component  velocity 
terms  have  been  normalized  by  M  . 


Table  A-2:  LDV  Velocity  Errors 


Ev- 

E«- 

Ev' 

8.1  % 

8.1  % 

2.0% 

8.0% 
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A.2  Total  Pressure  and  Temperature  Errors 

The  pressure  and  temperature  were  also  recorded  for  each  traverse  location.  The  error  in 
stagnation  pressure  was  estimated  as  0.3  %  based  on  the  fluctuation  from  the  mean.  Also,  the 
maximum  error  in  total  temperature  was  estimated  as  4.0  %.  This  estimate  was  based  on  the 
accuracy  of  the  measured  temperature,  and  the  assumption  of  constant  total  temperature  across 
the  boundary  layer. 

A.3  Calculated  Parameter  Errors 

With  the  errors  in  the  measured  quantities  (velocity,  temperature,  and  pressure) 
determined,  the  next  step  was  to  estimate  the  errors  of  the  calculated  quantities,  such  as  Mach 
number,  intermittency,  skewness,  and  boundary  layer  thickness.  By  substituting  the  previously 
defined  errors  into  the  data  reduction  equations  of  Chapter  6,  the  errors  for  the  calculated 
quantities  were  determined.  These  errors,  which  were  determined  by  the  error  analysis  of  Ref, 
[23],  are  summarized  in  Table  A-3. 
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Variable 

%  Error 

. 

u 

5.1 

V 

8.1 

u' 

2.0 

v' 

8.0 

u'v' 

8.2 

Q 

8.0 

M 

8.2 

Po 

0.6 

0.1 

z 

4.0 

28.0 

28.0 

Table  A-3:  Error  Summary 


X  17.0 

W 

dujdz  12.4 

dvidz  12.0 


Variable 

%  Error 

^RMS 

2.0 

RMS 

8.0 

SK 

8.5 

SK 

33.9 

Flu 

11.3 

FK 

45.3 

Y„ 

11.3 

Kv 

11.6 

duldx 

14.0 

dvidx 

14.0 

dujdy 

16.3 

dvjdy 

16.3 

Appendix  B:  Data  Files 


This  section  includes  a  listing  of  the  data  sets  from  the  present  study.  Included  are  the 
flat  plate  and  APG  processed  data  files  used  to  generate  the  figures  throughout  this  document. 
Table  B-1  presents  which  data  sets,  given  as  Tables  B-2  through  B-8,  were  used  in  each  figure. 
For  example,  the  data  files  listed  in  Tables  B-3  and  B-4  were  used  in  Figure  7-6.  Due  to  the 
number  of  measurement  stations  for  the  flat  plate  case  at  x  =  7 1  cm,  the  data  set  has  been 
presented  in  two  separate  tables  (Tables  B-3  and  B-4).  The  y/6„  data  column  has  been 
repeated  in  Table  B-4  for  convenient  comparison.  This  section  satisfies  the  criteria  of  Settles  and 
Dodson  [26],  by  providing  documented  and  tabulated  data  in  machine-readable  form. 


Table  B-1;  Tabulated  Data  by  Figure  Number 
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6.06E-05  4.4106  27.333  2522.6 


Appendix  C:  LDV  and  Flow  Field  Quality  Analysis 


C.l  8  vs.  16  Sample  Record  Length 

While  learning  to  use  the  LDV  system,  various  parameters  were  adjusted  to  determine 
which  BSA  settings  would  optimize  the  analysis.  One  of  these  adjustable,  and  user  defined, 
parameters  was  the  record  length.  Recall,  the  record  length  specified  the  number  of  samples 
used  to  generate  a  Fast-Fourier  Transform  of  the  Doppler  signal.  Each  BSA  was  restricted  to 
record  lengths  of  8,  16,  32,  and  64  [7].  Record  lengths  of  32  and  64  samples  were  not  used  in  the 
present  study  because  the  time  required  to  sample  the  signal  exceeded  the  traverse  time  of  the 
particle.  Thus,  because  of  the  high  velocities  the  analysis  was  restricted  to  either  8  or  16 
samples.  The  results  of  flat  plate  LDV  measurements  made  at  the  same  location  with  both  8  and 
16  sample  record  lengths  are  presented  in  Figures  C-1  and  C-2. 

Both  record  lengths  gave  close  agreement  on  the  mean  flow  measurements.  The  Mach 
number  profile  for  both  the  8  and  16  sample  settings  were  very  close,  with  the  8  sample  record 
length  predicting  slightly  lower  velocities,  and  hence,  slightly  lower  Mach  number.  The  largest 
difference  between  the  two  record  lengths  was  discovered  upon  comparison  of  the  turbulence 
parameters.  Figure  C-2  shows  a  comparison  of  the  u-component  turbulence  intensities  for  the 
flat  plate.  There  is  close  agreement  between  the  profiles  below  y/5  =  0.4.  As  the  measurement 
station  approached  the  freestream,  the  8  sample  case  predicted  freestream  turbulence  intensities 
that  were  three  times  those  given  by  the  16  sample  case.  This  difference  may  be  attributed  to 
noise  levels  that  were  too  high.  Based  on  this  analysis,  only  16  sample  record  length  LDV 
measurements  were  made  after  this  discovery.  The  8  sample  data  is  only  presented  in  this 
appendix. 
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C.2  Raw  vs.  Coincident  Data 


As  was  mentioned  in  Chapter  3,  coincidence  filtering  is  a  method  of  validating  the 
Doppler  signals.  Coincidence  filtering  determines  whether  or  not  the  BSA’s  are  tracking  the 
same  particle.  Raw  data  is  the  data  that  was  collected  by  the  BSA’s.  Coincident  data  is  that 
which  has  passed  the  coincidence  filtering  criteria.  Thus,  the  coincident  data  is  a  subset  of  the 
raw  data.  Therefore,  the  raw  and  coincident  data  were  compared  to  determine  if  the  coincident 
data  were  accurately  representing  the  raw  data. 

Figures  C-3  and  C-4  present  Mach  number  and  u-component  turbulence  intensity 
comparison  between  the  raw  and  coincident  data.  The  profiles  for  both  data  are  very  similar. 

From  the  raw  data,  all  turbulence  statistics  except  u'v'  are  calculated  by  the  BURSTware. 
Because  this  term  appears  in  the  turbulent  shear  stress,  it  was  necessary  that  the  coincident  data 
be  used  in  the  data  reduction  procedures  of  Chapter  6.  In  these  figures,  it  appears  as  if  the  raw 
data  LDV  measurements  were  made  closer  to  the  wall  than  the  coincident  data.  Actually, 
coincident  data  could  not  be  calculated  for  the  single  component  LDV  measurements  because 
only  one  component  of  velocity  was  measured,  and  the  coincidence  filtering  process  required 
comparison  of  the  Doppler  bursts  received  by  two  separate  BSA’s.  These  figures  show  that,  in 
fact,  the  coincident  data  have  the  same  trends  as  the  raw  data.  Therefore,  the  coincident  data 
may  be  used  to  present  the  processed  results. 


C-3 


1.6 

1.4 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

1.0000  1.5000  2.0000  2.5000  3.0000 

Mach  Number 


J - 1 - 1 - 1 - 1 - 1 _ I _ I _ ■  ■  I _ I _ I _ L. 


I  Coincident  Data,  APG,  68  cm 
I  I  Raw  Data,  APG,  68  cm 


n 


□i 


□ 


□  I 


D 


□ 


□ 


Figure  C-3:  Raw  vs.  Coincident  Mach  Number 
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Figure  C-4:  Raw  vs.  Coincident  Turbulence  Intensities 


C.3  Two-Dimensionality  check 

Another  criteria  for  Ref.  [26]  is  validation  of  two-dimensionality  of  the  flow  field.  If  the 
spanwise  flow  variations  ( duldz.,dvl dz )  are  smaller  than  the  order  of  the  other  strain  rates,  then 
two  dimensionality  of  the  flow  field  is  verified.  Strain  rates  in  the  spanwise,  or  w-direction,  were 
in  fact  smaller  than  in  the  other  directions. 

LDV  measurements  were  made  at  5  mm  on  either  side  of  the  tunnel  centerline  for  the 
strain  rate  measurements.  These  results  are  presented  in  Figures  C-5  and  C-6.  The  comparison 


of  the  fluctuating  cross-correlation  term,  u'v' ,  which  was  very  sensitive  to  variations  in  the 
flowfield,  is  very  good.  In  both  cases,  the  agreement  is  remarkable.  Thus,  the  flow  was 


considered  two-dimensional. 


Figure  C-5:  2D  Comparison  of  Velocity  Profiles 
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Figure  C-6:  2D  Comparison  of  Turbulent  Shear  Stresses 

C.  4  Repeatability 

Two  methods  were  used  to  determine  the  repeatability  of  the  present  study.  The  first 
method  was  to  compare  different  LDV  measurement  profiles,  which  were  measured  days,  and 
even  weeks  apart.  These  x  =  71  cm  flat  plate  profiles  are  presented  in  Figures  C-7  and  C-8.  The 

agreement  between  the  four  test  runs  is  very  good,  especially  for  the  —u'v'ju^  term,  which  is 
very  sensitive  to  flow  variations. 

For  the  second  method  of  repeatability  validation,  confidence  intervals  for  several  of  the 
mean  flow  and  turbulence  statistics  were  calculated.  Twenty  sets  of  data  were  taken  at  the  same 
traverse  location.  Two  of  these  sets  were  discarded  because  the  mean  velocity  was  much  lower 
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than  the  remaining  mean  velocities.  This  probably  occurred  because  data  acquisition  began 
before  the  tunnel  had  reached  steady  state. 


A  confidence  interval  is  a  measure  of  the  probability  that  the  error  in  the  measured  mean, 
X  ,  will  be  less  than  a  certain  value  when  compared  to  the  actual  value,  |J.j.  [2],  Thus,  a  95% 
confidence  interval  of  5%  indicates  that  there  is  a  95%  probability  that 


Because  the  value  |a,^ 
be  directly  calculated. 


<0.05 


(C.l) 


I  ^  ! 

was  not  known,  the  value  of  the  confidence  interval  in  Eq.  (C.l)  could  not 
Thus,  to  Student’s  t-distribution[2]  was  used: 


^  ^^n,a/2 

~  x^^N 


(C.2) 


where  N  was  the  number  of  data  points  to  calculate  x  ,  and  s  was  the  standard  deviation  of  x . 
The  t-distribution  was  a  tabulated  function  based  on  n  =  - 1  and  a  ,  which  described  the 

accuracy  of  the  confidence  interval.  Thus,  for  a  99%  confidence  interval,  a  =  0.01 .  The  results 
of  the  99%  confidence  interval  for  the  mean  flow  and  turbulence  statistics  of  the  18  valid  data 
sets  are  presented  in  Table  C-1. 


Table  C-1:  Repeatability  Confidence  Intervals 


Variable 

Confidence  Interval 

u 

0.15% 

V 

0.076% 

(u') 

\  /  rms 

1.53% 

■OH 

0.99% 

u'v' 

13.6% 

Variable 

Confidence  Interval 

u  Skewness 

2.60% 

u  Flatness 

3.23% 

V  Skewness 

14.1% 

V  Flatness 

8.00% 
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Appendix  D:  CFD  Validation 


As  mentioned  in  Chapter  1 ,  accurate  empirical  data  is  required  to  help  validate 
turbulence  models.  The  experimental  LDV  data  of  the  present  study  were  used  by  Fick[14]  to 
validate  a  compressible  turbulence  model.  He  numerically  studied  the  influence  of  pressure 
gradients  on  Mach  3.0  turbulent  boundary  layers  generated  by  the  same  geometry  as  in  the 
present  investigation.  Baldwin-Lomax  and  k-to  turbulence  models  were  used  in  a  cell-centered 
finite  volume  flow  solver  to  analyze  the  turbulent  flow  structure.  Agreement  between  the 
experimental  and  numerical  results  of  the  k-m  model  was  excellent  for  the  flat  plate  and  adverse 
pressure  gradient  cases[14].  The  numerical  results  are  presented  in  Figures  D-1  and  D-2.  The 
numerical  analysis  predicted  the  velocity  profile  very  well.  The  CFD  analysis  underpredicted  the 
turbulent  shear  stress  between  40%  and  80%  of  the  boundary  layer  thickness. 


Figure  D-1:  LDV  vs.  CFD  Velocity  Profile 
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